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Abstract. The use of the sine-Gordon equation as a model of magnetic flux propagation in Josephson 
junctions motivates studying the initial-value problem for this equation in the semiclassical limit in which 
the dispersion parameter e tends to zero. Assuming natural initial data having the profile of a moving —2tt 
kink at time zero, we analytically calculate the scattering data of this completely integrable Cauchy problem 
for all e > sufficiently small, and further we invert the scattering transform to calculate the solution for 
a sequence of arbitrarily small e. This sequence of exact solutions is analogous to that of the well-known 
A'^-soliton (or higher-order soliton) solutions of the focusing nonlinear Schrodinger equation. Plots of exact 
solutions for small e reveal certain features that emerge in the semiclassical limit. For example, in the limit 
e — > one observes the appearance of nonlinear caustics, i.e. curves in space-time that are independent of e 
but vary with the initial data and that separate regions in which the solution is expected to have different 
numbers of nonlinear phases. 

In the appendices we give a self contained account of the Cauchy problem from the perspectives of both 
inverse scattering and classical analysis (Picard iteration). Specifically, Appendix [A] contains a complete 
formulation of the inverse-scattering method for generic L^-Sobolev initial data, and Appendix [b| establishes 
the well-posedness for L?" -Sobolev initial data (which in particular completely justifies the inverse-scattering 
analysis in Appendix [A|. 



1. Introduction 

The sine-Gordon equation 
(1) e^utt - + sin(M) = 

describes a broad array of physical and mathematical phenomena. The partial differential equation ([T|) may 
be regarded as the continuum limit of a chain of pendula subject to an external (gravity) force and coupled 
to their nearest neighbors via Hooke's law. In nonlinear optics, the sine-Gordon equation is a special case of 
the Maxwell-Bloch equations and describes self-induced transparency in the sharp-line limit |16j . In biology, 
the sine-Gordon equation models transcription and denaturation in DNA molecules |18j . Backlund showed 
a correspondence between solutions of the sine-Gordon equation and surfaces of constant negative curvature 

i- 

In solid-state physics, the sine-Gordon equation models idealized magnetic flux propagation along 
the insulating barrier between two superconductors in a Josephson junction. Here the length of the 
transmission line corresponds to a length of 1 in terms of the dimensionless coordinate x measuring distance 
along the junction. Let L be the inductance per unit length and C be the capacitance per unit length. 
Then v := (LC)"^/^ is the typical velocity parameter, and the macroscopic time scale t measures one 
dimensionless unit when t^/v seconds have passed. The parameter e is the ratio of the Josephson length 

1/2 

to the transmission line length ^o- The Josephson length ij is in turn proportional to <i>Q , where 
<l>o := h/ (2e) « 2.064 x lO"""^^ V sec is the quantum unit of magnetic flux. Laboratory experiments by Scott, 
Chu, and Reible |20| analyzed flux propagation in Josephson junctions of length = 35 cm for which ij 
was approximately lO""* to 10"'^ m. Therefore, in these experiments, e := Ij/Iq ~ 0.0005. The period of a 
signal input to the transmission line in these experiments was typically on the order of 10~^ seconds, which 
is approximately one dimensionless time unit on the i-scale. Together with e being small, this motivates 
the study of the semiclassical (or zero-dispersion) limit as e | 0. For analytical convenience we choose to 
study the Cauchy initial-value problem on the real line a: € M. Formulating a semiclassical Cauchy problem 
means fixing functions / and g independent of e, and then, for all e > sufficiently small, posing the Cauchy 
problem for (jlj with initial data of the form m(x, 0; e) — f{x), eut{x,0; e) — g{x). See Appendix [b] for an 
account of the well-posedness theory of this Cauchy problem for e > fixed. Solving the semiclassical Cauchy 
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problem means obtaining the one-parameter family of solutions u(x, t; e). We are usually most interested in 
the asymptotic behavior of the solution as e | 0. 

In this paper, we consider the sine-Gordon equation (jlj for all e > sufficiently small with the 
initial condition 

(2) ^™ (^2) •==sech(a::), := tanh(a;), g :— 2fisech{x), x G M. 

where /i G M is a parameter. We refer to the solution of the Cauchy problem as u{x , t; e , fi) . The topological 
charge (or winding number) of solutions satisfying ([2| is a constant of motion given by 

1 f+°° 

(3) := — / Uxdx = -1. 

^TT J-00 

From one point of view, the initial data Q are natural to study, because u{x,t; ^^^^ + 1, /i) is an exact 
mathematical antikink solution of the sine-Gordon equation explicitly given by 



(4) 



i{u(x,t; \/^'^ + 1,^)) = 1 — 2sech^ ( — {x — vt)] , 

\ew ) 

i{u(x,t; \/ ij'^ + 1, 11)) — 2 sech [ — (x — vt) ) tanh [ — (x — vt)\ 

\£w J \ew J 



with velocity v and width parameter w given by 

(5) — ; ^ / „ - 1. w = -\^/l-v'^. 

From these formulae we see that u is a traveling wave with velocity v bounded by 1 (the light speed), 
demonstrating the hyperbolicity of the sine-Gordon equation. This solution admits a natural relativistic 
interpretation since the relationship between v and w corresponds to Lorentz contraction in special relativity. 

For £ 7^ fi^ + I, the initial data Q no longer corresponds to simply one soliton, but in general 
excites a nonlinear superposition of kinks, antikinks, breathers, and radiation. It is interesting to observe 
that the initial data ([2| satisfies the advection equation ut + vu^ = with constant velocity v = fx/e. In 
this sense, we may consider the initial data as being in uniform motion to the right with velocity v. Note, 
however, that if /i 0, then for e > sufficiently small, the velocity v of the initial data exceeds the 
constraint \v\ < 1 imposed by the hyperbolic nature of the sine-Gordon equation ([l]). In this situation, one 
might expect the sine-Gordon equation to regularize the superluminal velocity of the initial data for t > by 
some kind of catastrophic effect that destroys the profile of the initial data. In fact, we will show (see figures 
[6] and [S]) that the regularization of the velocity takes place via the emission of a large number (inversely 
proportional to e) of kink-antikink pairs. 

The family of solutions corresponding to the initial data ([2]) may be viewed as an analogue for 
the sine-Gordon equation of the A^-soliton (or higher-order soliton) solution to the cubic focusing nonlinear 
Schrodinger (NLS) equation 

(6) iqt + ]^qxx + \q\^q^Q ■ 

Satsuma and Yajima [H] found that with initial data q{x,Q]N)—N sech{x) the scattering data relevant for 
the focusing NLS equation can be found in closed form for any G M. Furthermore, if G Z then the 
scattering data are reflectionless and so the solution q{x, t; N) can be found more-or-less explicitly. In |T7j it 
was noted that, with e = c/N and r = Nt/c, the function (p{x, r; e) = cq{x, t; N)/N satisfies the initial- value 
problem 

(7) ie(t)r+—<j)xx + \4'\'^(f>^0, (l){x,0;e) ^ csech{x). 

The functions 4>{x,t;£) therefore solve a semiclassical Cauchy problem since (/>(x, 0; e) is independent of e. 
Numerical reconstruction of the inverse-scattering solution for e = En = c/N, N — 1,2,3, . . . in [17] revealed 
a spatio-temporal pattern for (j) emerging as e | consisting of a fixed macrostructure with nonlinear caustics 
(phase transition boundaries or breaking curves) separating regions of the space-time plane consisting of 
oscillations of different local genus (number of nonlinear phases). At least two caustic curves appear in 
the dynamics (a primary caustic t — ti(x) and a secondary caustic t ~ ^2(2^) > ti(x)). The semiclassical 
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asymptotics for times t up to and just beyond the primary caustic were obtained in [12J and these results 
were extended to times t just beyond the secondary caustic (requiring a substantial modification of the 
method that captures the primary caustic) in [T^. In a related result, Tovbis and Venakides ^2T] generalized 
the calculation of Satsuma and Yajima by computing the scattering data associated with the semiclassically- 
scaled focusing NLS equation Q explicitly for all sufficiently small e > when the initial data is given in 
the form 

(8) 0; n) = A(a;)e'■5(^)/^ S'{x) ^ -^tanh(x), A = csech(a;) . 

Subsequently, the Cauchy problem for ([t]) with this initial data has been studied by Tovbis, Venakides, and 
Zhou In this paper, we present a calculation of the scattering data for ([T|-([2]) analogous to the 

work in [19] and [21], and we also present an explicit computation of u{x,t;e,fi) as e | analogous to [TTI. 
The asymptotic analysis of the semiclassical Cauchy problem for sine-Gordon corresponding to the work in 
[m [121 [551 [23] will be carried out in a later work. 

The sine-Gordon equation ([T]) is an integrable system, possessing a Lax pair (see (123) and (124l) 



and admitting all the benefits thereof, including the existence of inverse-scattering transforms for solving 
Cauchy problems in various coordinate systems. We consider the Cauchy problem in laboratory coordinates 
and we use the Riemann-Hilbert formulation of inverse scattering. For the sine-Gordon equation in char- 
acteristic coordinates, the inverse-scattering method was first given in [T] and [2S]. The inverse-scattering 
method corresponding to the (noncharacteristic) Cauchy problem for the sine-Gordon equation in laboratory 
coordinates was worked out by Kaup [13j . An account of the Riemann-Hilbert method for carrying out the 
inverse step in laboratory coordinates can be found in the text of Faddeev and Takhtajan [8], and further 
developments to the theory were made by Zhou [26] and Cheng et al. In our paper, we add to this lit- 
erature by giving in Appendix [X] a complete description of the Riemann-Hilbert formulation of the solution 
of the Cauchy problem in laboratory coordinates assuming that at each instant of time the solution u has 
L^-Sobolev regularity. That the sine-Gordon equation ([TJ preserves this degree of regularity if it is present 
at i = is established by independent arguments in Appendix [B] 

Briefiy, the inverse-scattering method proceeds as follows. Cauchy data for the sine-Gordon equation 
characterize a set of scattering data, which consist of the reflection coefficient p : K ^ C, the eigenvalues 
{zn}, and the modified proportionality constants {c„}. The scattering data are used to formulate a Riemann- 
Hilbert problem with an explicit, elementary dependence on x and t. While it is not in general possible to 
solve a Riemann-Hilbert problem in closed form, for reflectionless Cauchy data (i.e. for which p(z) = 0) the 
Riemann-Hilbert problem can be reduced to the solution of a system of linear algebraic equations. 

In Section [2] we explicitly calculate the scattering data corresponding to viewing the initial data 
([2]) as a kind of potential in the linear scattering problem ( 123 1 associated with the Cauchy problem for the 



sine-Gordon equation ([ij. Our analysis will be valid for all /j, G K and e > sufficiently small. Furthermore, 
we show that if e lies in the sequence 



(9) e = eN{fi):=^^^, iV e Z+ {0, 1, 2, . . . } 

(note that this sequence converges to zero as ^ 00), then the scattering data are reflectionless ultimately 
implying via inverse-scattering theory that u can be computed explicitly (that is, u can be expressed by a 
finite number of arithmetic operations). The inverse step is carried out for e in the sequence ^ corresponding 
to reflectionless initial data in Section[3] where cos(u) and sin(M) are extracted by considering an appropriate 
limit of the solution to the Riemann-Hilbert problem. As e | through this sequence, a pattern emerges in 
which u consists of modulated wave trains of wave number and frequency inversely proportional to e with 
one or more nonlinear phases. The spatio-temporal scale of the modulation is fixed as e | 0. Regions of 
space-time containing waves with different numbers of nonlinear phases are separated by nonlinear caustics 
that are independent of e for fixed fi. See figures [5] |6] and [s] for plots of cos(u) exhibiting these features for 
various values of e and /i. At a qualitative level, these features resemble those observed for solutions of the 
semiclassical Cauchy problem for the focusing NLS equation. Section [4] discusses the limitations inherent in 
an approach to the semiclassical limit based upon calculations of complexity and sensitivity increasing with 
N ^ e^^, and explores possible extensions. 
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Remark 1. Under the scalings x — eX and t — eT and the choice e — eN{^) (see ([9|), equations Q and ^ 
become 



(10) 



U' 



TT 



u 



XX 



sm{U) = 0, 



(11) 



-UiX,0;N,^i) ] =sech(ejv(MW 



cos 



-U{X, 0; N,n)]= tanh ieN{^i)X) , 



C/T(XO;iV,^) = 2^sech(ejv(M)^) , 
where U{X,T; N, fi) — u{x , t; e , fi) . This is a fixed-dispersion Cauchy problem with a sequence of different 
initial conditions depending on N, just as in the problem for the NLS equation studied by Satsuma and 
Yajima. The initial conditions all have topological charge w[U] = —1 but U{X,0; N, fi) and Ut{X,0; N, ii) 
become more dilated in X (slowly-varying) as N increases. Therefore, an alternate way of viewing our result 
is that we can find exact solutions to the fixed-dispersion initial- value problem ( 10 1-( 11 ) for N G Z+. As an 
example of an explicit solution of ( |10[ ) obtained in this way, when N = 1 and /i 

2n{X, Tf 



we have 



(12) 
where 

(13) 



n{X,T) 



11 



cos 



d{X,T) 4 



2 + 3 cos 




cos(C/(X, T- 1,0)) 



cosh 



cosh 




2 cosh 



^-X 
3 



cosh 



\x 



cosh(X) + 2 cosh -X 



The focusing NLS equation (|6| admits a scaling symmetry in which scaling the independent variable x is 
equivalent to scaling the dependent variable (amplitude) q and the time t. Thus, the iV-soliton (or higher- 
order soliton) solutions of the focusing NLS equation that were originally obtained by Satsuma and Yajima 
[1.9j by considering a fixed-width pulse with variable amplitude can just as easily be viewed as a fixed- 
amplitude pulse with variable width. From the point of view of semiclassical asymptotics, dilation in x is the 
more natural interpretation of the higher-order solitons as the presence of the parameter e in ([T]) amounts to 
rescaling x and t, and thus the variable width of the pulse is absorbed into the semiclassical parameter e as 
in |17| |12j . Of course, the sine-Gordon equation does not admit the amplitude/dilation symmetry enjoyed 
by the focusing NLS equation, so we are not free to interpret the family of exact solutions we obtain in this 
paper in terms of scaling of amplitude. It seems that perhaps a more generally fruitful approach to seeking 
analogues of the higher-order soliton in other integrable systems is to consider pulse width dilation as being 
more fundamental than amplitude dilation. As more evidence of the utility of this approach (beyond the 
sine-Gordon example), the modified NLS equation (which includes an additional term in ^ that breaks the 
scaling symmetry needed to exchange amplitude for width) does not have higher-order solitons in the sense 
of Satsuma and Yajima, but it does have exact solutions corresponding to arbitrarily width-dilated pulses 
that are useful in semiclassical analysis [5]. 



Remark 2. In characteristic or light-cone coordinates x and r defined by x = x -I- t and t = X — t, the 
sine-Gordon equation (jlj is e^v^r — sin(w), where w(x, e) — u{x, t; e). The associated x evolution equation 
in the Lax pair is the Zakharov-Shabat eigenvalue equation, which is the same eigenvalue equation as for the 
focusing NLS equation [24 . Thus it is possible to solve a semiclassical characteristic Cauchy problem with 
special initial data w(x,0;e) — csech(x) using the Satsuma- Yajima higher-order soliton solution. However, 
in many applications (as in Josephson junction theory) , the correct problem is the non-characteristic Cauchy 
problem with two independent initial conditions: m(x, 0; e) = f{x), eut{x, 0; e) = g{x). The Satsuma- Yajima 
solution to the semiclassical problem posed along a characteristic t = Q ov x = t will have a very complicated 
and unwieldy form and an undesired dependence on e upon restriction to t = 0, and therefore is probably 
not relevant to the non-characteristic semiclassical Cauchy problem we wish to consider. 

On notation. As will be explained in detail in Section [2] and Appendix [A] we will use three different gauges 
for the eigenvalue problem. Objects associated with the infinity gauge will be denoted by an overline (J). 
Objects associated with the zero gauge will be denoted by an underline (J). Finally, objects associated with 
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the symmetric gauge will not have a bar (J). The complex conjugate of a is denoted by a*. We use the 
notation /(a, 6; c) to emphasize that c is a parameter. The dependence on parameters may be suppressed by 
writing /(a, b) in place of /(a, b; c). We also make frequent use of the standard Pauli matrices defined as 



(14) 



0-1 := 



"0 


1' 









'1 " 


1 





, era := 


i 


, 0-3 := 


-1 



Vectors will be denoted by bold lower-case letters and matrices by bold upper-case letters, with the exception 
of the Pauli matrices. The transpose of a vector v is denoted by v^, and the conjugate-transpose of a matrix 
A is denoted by A^^. Finally, xs{') denotes the characteristic function (indicator function) of a set S, that 
is xs{x) = 1 if .T e S" and xs{x) ~ otherwise. 

2. Scattering Theory for the Special Initial Data 

The quantities z zt appear throughout the scattering and inverse-scattering theory of the sine- 
Gordon equation ([T|), and so for convenience we define 

(15) D{z):=l(^z+^y E{z):^\(^z~^ 

Fix the initial condition Q. To find the scattering data it is necessary to solve the following generalized 
eigenvalue problem for 3{x; z,t,e,^) (see, for example, [TS]): 

4:E + z~^{l — cos(w)) —z~^ sin(u) — ie{uj. + Ut) 
-z^^ sin(M) + i£{ux -t- Ut) —AE — z^^{\ — cos(w)) 



(16) 



This formulation of the eigenvalue problem is useful in the study of solutions J when z is bounded away 



from (see [13] and Proposition A.l below), and for this reason, we say that (16 1 is written in the infinity 
gauge. The use of alternate gauges proves to be beneficial. For example, the gauge transformation (1391 
(see Appendix [A| casts the eigenvalue problem (16 1 into an alternate form that is useful in the analysis of 
solutions corresponding to bounded z, and in particular near z = (see [13] and Proposition A.2|. Therefore, 
we refer to the coordinate system arrived at via the transformation ( [139| ) as the zero gauge. While the infinity 
gauge and the zero gauge are useful in the analysis of the scattering problem required to formulate an inverse- 
scattering theory, to calculate the scattering data corresponding to (pi we found it to be useful to introduce 
a gauge transformation that symmetrizes the appearance of z and z in the eigenvalue problem and at the 
same time also removes the function Ux from the coefficients. It is in this third, symmetric gauge that it is 
easiest to see the eigenvalue problem is in fact hypergeometric for the initial data ([2]). Once it is clear from 
working in the symmetric gauge that the eigenvalue problem has exactly three regular singular points, it is 
possible to use the theory of Euler transforms to analyze the asymptotic behavior of the Jost solutions and 
thus obtain the scattering data. 



2.1. Transformation to a hypergeometric equation. The first step in transforming (16 1 into a hyper- 
geometric equation is to introduce an appropriate gauge transformation. If J satisfies equation (16), then 
the invertible transformation (having an interpretation as a rotation at each x by an angle —u/A) 



(17) 



J(a;; z, t, e, ^) = AJ 



yields a solution J of the eigenvalue problem written in the symmetric gauge: 

ieut 



(18) 



AieJx 



AEcos f " 
V2 



-4D 



-AD sin ( - 
^ 2 



-4£:cos - 



(i) 



Written in this forn:Q the eigenvalue problem appears similar to one used by Faddeev and Takhtajan (see 
[S] part 2, chapter 2, equation 4.1). The Jost solutions J* in the symmetric gauge are defined to be the 

^The absence of Ux in the symmetrized form l |18| of the eigenvalue problem provides an alternate framework in which to 
consider discontinuous initial data without the use of delta functions (c/. ITT\ ). 



fundamental solution matrices of the linear problem (18 1 for real values of z, normalized by the conditions 



(19) 



J+ = 





^(,-iEx/e 
'^-iEx/e Q 



+ o(l) as X 



'Oo for z e 







piExje 



+ o(l) as a; ^ +oo for z G 



We denote the columns in this way: J* =: [jf They are related to the Jost solutions for the infinity 
gauge (see (125l) by 

(20) jj;=Aj2", j2" = Aj7, j+ = Aj^, j+ = Aj^. 

For the choice of initial condition (|2]) , the eigenvalue equation ([T8| takes the form 

2£;tanh(x) (-21) - i/z) sech(a;)" 

(-2D + i/i) sech(a;) -2E tanh(a;) 



(21) 



2jeJT. = 



J. 



With the change of independent variable 
(22) 



y = tanh(a;), 



the eigenvalue problem (21) becomes 
(23) 2ze(l - 



2Ey {-2D - ifi) {I - y'^y^'^' 

(-2L> + i/x)(l-j/2)i/2 ^2Ey 



Here — 1 < y < 1 and the positive square root is chosen. There are two ways to eliminate the square roots 
in the coefficient matrix. The first is to introduce the linear transformation 



(24) 



Ji 



1 

(1-2/2)1/2. 



which results in a differential equation satisfied by g{y; z, e, /i): 

2Ey {-2D - i/i) (1 - 2/2) 
-2D + iyL {-2E + 2ie)y 



(25) 2*£(1 - 2/2)gy = 

This equation has exactly three regular singular points 2/G{— l,l,oo} and can be written in hypergeometric 
form. We will use (25 1 to find expressions for and jj^. 



Remark. If we had taken (16 1 instead of (18 1 as our starting point and followed analogous steps, namely 
(i) substitution of the initial data using double-angle formulae, (ii) the independent variable transformation 
y = tanh(a;), and (iii) the use of the gauge transformation (24 1 to reduce the problem to rational form, we 
would have arrived at 

'2E + z-\l - 2/2) -z-i(l - j/2)y + ie{l - y'^y 
-z^^y-ie -2E - z^'^{l - y'^) +2iey 



(26) 



2ie{l - 2/2)gy 



instead of (251. Let v = y~^. Then near v = 0, (26) has the leading-order form — 0{v~'^) ■ g, whereas 
(25 1 has the leading-order form gy — 0{v^^) ■ g. After some calculation it is possible to see that the method 
of Frobenius still applies to ( 26 1 near v — even with the additional growth ai y = oo due to special 
identities holding among the entries of the matrix coefficients of the leading-order terms. However, the local 
(Frobenius) analysis of (25) is more straightforward with only a double pole at u = 0. Later we will also see 
that it is more difficult to obtain integral representations for solutions of ( 26 1 than for ( 25 1 . 



(27) 



An alternative to the linear transformation (|24l is 





1 



(1-2/2) 



-1/2 



which after substitution into (23) yields the differential equation for h{y; z,e,fi): 
(28) 2ie{l-y^)hy 



2Ey -2D - ifi 

{-2D + {I - y^) {-2E-2ie)y 
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This equation also has exactly three regular singular points j/g{ — l,l,oo}. It will be used to find expressions 
for j+ and j2". 

2.2. Integral representations for Jost solutions. We now use the theory of Euler transforms [TU] to 
derive integral representations for the four Jost solutions, starting with and j^. Define 



(29) 7 = 7(m) Va^' + I- 

Proposition 2.1. Choose the principal branches of the functions (s ± l)~i/2+7/2e mith branch cuts on the 
real s-axis from ^1 to — oo. Also choose the principal branch of {s — y)"'^/^~T/2'^~i with branch cut on the 
real s-axis from y to —oo. Take to be a closed contour in the s-plane passing through the branch point 
s = —1 and encircling s = 1 once in the counterclockwise direction (see figure 1(a)). Then, for z G M\{0}, 





(a) The contour S 



(b) The contour S" 



Figure 1. The integration contours and branch cuts for the Jost functions. 
)n ^- [Jii^JtiV «s by 

(30) J+ =Ci(l-y2).B/2e I (s_l)-l/2+7/2e(^^^)-l/2+W2e(^_y)-.i=;/e-7/2s^^ 

(31) J2+ = Ci(l - y2y£/2.+ l/2 A-(7 + /^)A /■ _ ^)-l/2+7/2e(^ ^ ^^_l/2+V2e(^ _ yy^E / e--, /2e-l 

+ (7 + WV JE+ 

with 



(32) 



C'l = Ci(z,e,/i) 



r(l + f)r(i + ^-f) 

2l+7/2e,,r(i + ^) 



Proof. We begin by computing the Frobenius exponents of ( p5| . Assume that, for some number p+ g C, g 
has a Frobenius series about y = 1 of the form 



(33) 



g(2/) = Eg»(i-y)'^ 



+« 



n=0 



for some vector- valued coefficients g^ . Substituting this series into equation ( 25 ) and considering the leading- 
order terms immediately gives the (indicial) eigenvalue equation 

2E 1 + 

-2D + in -2E + 2ie\ ^° ' 

Therefore, the Frobenius exponents at ?; = 1 are 

iE iE 1 



(34) 



4iep+g(j 



2e ' 2e 2 ' 



(35) 

Similarly, substituting a series of the form 

00 

(36) g{y) = Y.Sna + yr 

n=0 



+n 



into equation ( 25 ) and considering the leading-order terms shows that the Frobenius exponents ai y — —1 

iE 



are exactly the same: 
(37) 



P = 



iE 1 

2£' ^2e^2' 
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We now shift two of the exponents to zero via the substitution 
(38) g(y) = (l-y2)^^^/2ef(y). 

It follows that f(y; z,e,fj.) satisfies the differential equation 

{-2D - ^/i) (1 - y^) 



(39) 



2*6(1 - y'){y 



-2D + ifi {-AE + 2ie)y 



We attempt to express f —: [/i,/2]^ as 



(40) fi{y)^ F,{s){s-yrds, f2{y)= F2{s){s-y)^ds 

JT.+ JT.+ 

where the Euler transforms i^i,2(s; z, e, /i) and the constant exponents a — a{z,e,fj,) and (3 — /?(z,e,/i) 
remain to be chosen. Substituting the expressions (40 1 into the system (39) shows that the first equation 
2iefiy = {—2D — i^) /2 can be easily solved by choosing 

(41) ^^(^) = (2iJT7^^^(^) ^ = "-1- 

Remark. The fact that there is such a simple relationship between Fi and F2 is related to the fact that /2 
could be easily eliminated to write a second-order differential equation for /i that is essentially the Gauss 
hypergeometric equation. On the other hand, if we had worked from the beginning in the infinity gauge, 



the elimination of /2 using the first row of ( 26 ) would have resulted in a second-order equation that is not 



obviously of hypergeometric form. This, in turn, would lead to further complications in the following analysis 



leading from ( 43 1 to ( 49 ) 



(42) 



It remains to satisfy the second equation of the system ( 39 ) 

f2y 



2ie{l - y^)f2y = {-2D + z/i) A + {-AE + 2ie) y/a. 



Writing y = s — {s — y) and 1 — y^ — {1 — s^) -|- 2s{s — y) — {s — y)"^ and using equations (41 1 gives 

Fi(s) [(1 - s^){s - yr-^ + 2s{s - yr-^ - {s - y)"] ds 



-2ie'^a{a - 1) 



{2D + in) 



(43) 



^D+^^ 



2^ ./s+ 



F^{s){s-yrds 
2iea 



{2D + ifi)Js+ 



If we now choose a to satisfy the quadratic equation 

2i£'^a{a - 1) 



(44) 



{2D + i^i) 



^ iD 



- - {2iE + e) 

2 J ^ ' {2D + ifi) 



then the {s — y)" terms will cancel in equation (43). Specifically, we choose 

Using integration by parts to eliminate the (s — y)" ^ term yields 
d 



(46) 



4. d-s 



)] (s - y)"^^ ds + 2£{a - 1) / Fi{s)s{s - y)""^ ds 



E+ 



= {-2iE-e)[ Fi{s)s{s-yY 



Setting the integrands equal and using equation ( 45 1 for a gives the first-order linear differential equation 
(47) e{l-s^)^^F,{s)^{e--i)sF,{s) 



(48) 



ds 



logFi(s) 



1_ 
2e 



2 I s 
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1_ 

2e 



for Fi{s). The general solution 

(49) Fi{s) = Ci(s - l)-i/2+^/2^(,s + l)-i/2+7/2. ^ 



where Ci is an integration constant, gives equations (30 1 and (31) up to the choice of the constant Ci. 

T he constant Ci is to be chosen so that is normalized as required in equation (19 1. Consider 
equation (30) for Jj^ as y ^ 1. Now for s £ S+, 



(50) (S - l)-l/2+7/2.(5 i)-l/2+^/2e(^ _ yy,E/e-'t/2e 



< 



is- I) 



-l/2+7/2e 



-iE/e-j/2e 



-l/2+7/2e 



<i^|s+irl/2+7/2e 



for some constant K > 0, as is bounded away from s = y and s = 1. Since "f/2e > 0, the function 
\s + l|^i/2+7/2e jg integrable on E+. Therefore, by Lebesgue's dominated convergence theorem. 



(51) 



J+=C,il-yy^/'^ 



{S + l)-l/2+''/2-(s - l)--E/e-l/2 ^ 



as y ^ 1. 



The integrand is integrable at s = 1, so we can collapse S+ to the upper and lower edges of the branch cut 
["1, 1], yielding 



(52) 



(s+l)-l/2+7/2e(^_l)-.B/.-l/2^^ 



2zsin (^TT + *f ) ) /' (S + l)-l/2+7/2e(^ _ ,yl/2-^E/e^^ 

(s+l)'l/2+^/2e(^_^)-l/2-^£/e^^_ 



_ 2i7r 

"r(i + f)r(i-f) 

In the last step we used the reflection identity 

(53) sin(7rz) = ^ r. 

r(z)r(i - z) 

The remaining integral is a beta integral, which may be expressed in terms of gamma functions. Indeed, 
making the change of variables s = 2w — 1 gives 

{S + l)-l/2+7/2.(l _ ,yl/2-^E/e ^ 2^/2e-.£/s f' ^-l/2+7/2s(^ _ ^yl/2-^E/e 
L Jo 



(54) 

using the identity 
(55) 



_ 2l/2e-iE/i 



2e e ) 



£(i±l)rG-f) 
r(i + 

r(a)r(6) 

T{a + b) 



valid for ^{a),3l{h) > 0. Also note that, as y ^ 1, 



(56) 

Therefore, as a; — > +oo, 
(57) 



1 - y2 ^ sech'(a;) = 22e-2-(l + o(l)). 



r(i 



7 iE_^ 

2e e ) 



0(1) 



Comparing equations (19) and (57) gives the expression (32) for the constant Ci 



□ 



Proposition 2.2. Choose the principal branches of the functions (±1 — s)^i/2+7/2£ ujitk branch cuts on the 
real s-axis from ±1 to +oo. Also choose the principal branch of [y — s)^'^/*^ r/2£-i y^aii branch cut on the 
real s-axis from y to +oo. Take to be a closed contour in the s-plane passing through the branch point 
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s = 1 and encircling s = — 1 once in the counterclockwise direction (see figure 1(b)). Then, for z £ M\{0}, 
jr =: [-^n' J21V is given by 

(58) J-^ = Ci(l-y2)»i^/2e I (_i_^)-l/2+7/2s(i_^,)-l/2+-y/2e(y_^)-.£/e-7/2e^^ 

(59) 

V ^ + (7 + m)* 



(_1 _ s)-l/2+^/2-(l - s)-l/2+^/2s(y _ ^^y^E/ e-^/2e-l 



with Ci given by equation (32). 

Proof. The construction follows that of j^*", except with in place of E+, and choice of the solution 

(60) Fi(s) = Ci(-1 - .s)-l/2+7/2.(l _ 5)-l/2+7/2e 



to equation ( 48 1 . 



□ 



Now we use the system of differential equations in the form (28) resulting from the transformation (27 1 to 
find expressions for the Jost solutions and j^. 



Proposition 2.3. Take (s ± i)-i/2+7/2e^ (^ _ y)-i-E/s-7/2e-i^ j]+ Proposition T/ien, /or 

z e M\{0},. j+ =: [J+ , 73+2]"^ IS given by 



(61) J+ =-C2(l-2/2)-«B/2.+l/2 + y (g_l)-l/2+7/2e(^_^i)-l/2+7/2e(^_y)^i;/e-7/2e-l 



(62) 
with 
(63) 



J+ = C2(l - J/2)-^^=/2e / _ i)-l/2+7/2.(^ ^ ^)-l/2+7/2e(^ _ yyE/e-^/2e 



C2 = C2(z,e,/i) := 



r(^-f)r(i+2^ + f) 

2l+7/2e,,r (i + ^) 



Proof. The Frobenius exponents p around y = —\ and p+ around ?/ = 1 for h satisfying (28 1 are 

iE iE I , iE iE I 

(64) = — , h -, p+ = — , h -. 

^ ^ 2e ' 2e 2 ' ^ 2e ' 2e 2 

Therefore, defining f (y; z, e, /x) = [/i(y; z, e, /i), /2(?/; e, /i)]^ in terms of h by 

(65) h(y) = (l-y2)-s/2.+i/2f(j^) 

has the effect of shifting one exponent to zero near each of the points y = ±1. By direct calculation, f 
satisfies 



(66) 



2ie{l - 2/2)f^ 



{AE + 2i£) y -2D - i^jL 
{-2D + iy){\-y^) 



Assume integral representations of the form 



(67) 



/i(y)-/ Fi(s)(s~2/)"ds, f2{y)= F2{s){s~yYd 



Proceeding as in Proposition |2.1[ we obtain 

F,{S) ^ - ('l±Sri±}^\ F^(^,y FM = C-,{S ~ l)-V2+7/2e(, ^ i)-l/2+7/2 

\z- (j + fJ-)iJ 



(68) 



a = (3-l, P 



iE 7 



e 2e 

where C2 is a constant of integration. 

The constant C2 is chosen so j J is normalized as required in equation (19 1. Starting with equation 

( 62 1 for J^^ , for s £ we have 



(69) 



(.-1) 



-l/2+7/2e 



(^+1) 



-1/2+7/26^^ _ yyE/e-j/2e 
10 



<i^|s+ 11-1/2+7/25 



for some constant K > 0. By dominated convergence, 



(70) 



(S - + ly^/'-'/^ds + 0(1) 



as y 1 . 



This is the same expression as equation (51 ) for J^^ with Ci and E replaced with C2 and —E, respectively. 
Therefore, as y ^ 1, 



(71) 



^22 — <-"2e 



2l+7/2£j 



1 ^£;\\r(i + ^)r(i 



^ V2 e 



)r(i. 



2£ 



r(i + 5^ + f) 

o(l) 



0(1) 



using the identity ( 53 ) , which gives equation (|63 



□ 



Proposition 2.4. Z)e/?ne (±1 - s)-i/2+7/2£^ (^y _ g-^-tE/e-^/2e-i ^ andY.- as m ProposiUon\2^ Then, for 

z g M\{0}, ^2 ^- [Ji2^ 'J22V ^•s given by 
(72) 

,(1 - y2)-^£/2e+l/2 + y (_1 _ 5) - l/2+7/2e ( ^ _ ^) -l/2+7/2s _ ^).i;/e-7/2e-l ^ 



■"12 



-Co 



(73) 



^2(1 - y2)-»£/2s / _ ^)-l/2+7/2s(^ _ ^)-l/2+7/2e(y _ ^yE/e-j/2e 



with C2 given by equation (63). 

Proof. The integral representation for is derived in the same way as the representation for j J in Proposition 

□ 



2.3 



with S- replacing S+ and the choice J^2(s) = C'2(-l - s)-^/'^+^/'^%\ - s)-i/2+'^/2 



2.3. The scattering data. We now use the integral formulae from Section [273| to calculate the scattering 
matrix S, the eigenvalues and (in certain special cases of interest) the proportionality constants {rjn}. 

Proposition 2.5. The coefficient 6*22(2) — 322(2; e, fi) is given by 



(74) 



822(2) 



(^-(7 + M)^:) [rik~f)f 
(z + (7 + /i)*) r(i-2|-f)r( 



jr ii 

2e e 



z e 



Proof Take z e M\{0}. By equations ([T67| and ([20]), 

(75) = 622jr + ^12j2"- 

To determine 622 we use 

(76) ^22 = fl = - lim det = lim (j+e-'^^/^ 
We now analyze Jj^ as a; — > — 00. Consider the integral 

(77) h(y; Z, e, m) := / (s - l)-l/2+7/2e(^ + ^)-l/2+7/2e(^ _ yyE/e-y/2e-l 

JT.+ 

Making the substitution s — --X + (y + \)v gives 



(78) 



h = (y+ 1)^^/^-1/2 / ^ ^ ^^^^_i/2+^/2e^-l/2+7/2e(^ _ ^^.B/e-7/2e-l 



S + 



=: (j/ + 1)^^/^-1/2/2. 



Here we take S+ to be the counterclockwise-oriented contour starting at w = 0, following the semicircle in 
the lower half-plane of unit radius centered at w = 1 to w = 2, proceeding along the real axis to w = 2/(?/+ 1), 
coming back along the real axis to w = 2 along the top side of the branch cut, and t hen re turning to u = 
along the semicircle in the upper half-plane of unit radius centered at u = 1. See figure 2(a) Now for v G S+, 



(79) 



\-2+(y+l)v\ = \y+l\ 



y+l 



\y + M 
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s = —1 — i 

(a) The contour S+. (b) The contours Sl+ and E2+. 



Figure 2. The integration contours used to calculate 5*22(2) and Si2{z). 



The inequality above follows because the contour E+ lies inside the circle of radius 2/{y + 1) centered at 
V = 2/{y + 1). Therefore, for u G S+ with |i; - 1| = 1, 



(80) 



(_2 +(y+ l)«)-l/2+7/2e„-l/2+7/2s(^ _ ^yE/e-j/2e~l < ^|^|^l/2+7/2e 



for some constant K > 0. The right-hand side is integrable on the part of S+ with \v — 1| — 1. For the 
part of I2 involving integration over the upper and lower edges of the branch cut [2,2/(y + 1)], we have 
I - 2 + (y + l)v\ < \2y\ < 2 and X[2,2/(y+i)] («) < X[2,+oo)(w)- Therefore, since -1 + 7/2 > for e > 
sufficiently small, 

(81) \-2+{y+l)v\-'/'+''/'^X[2,2/iy+i)]{v) < 2-'/'+^/'^X[2.+oo]{v). 
So, forz;e [2,2/(y+l)], 

(82) (-2 + {y + l)i;)-l/2+7/2e^-l/2+7/2e(^ _ ^yE/e-^/2e-l < 2-l/2+7/2s^-l/2+7/2e(^, _ l)-^/2^-l, 

which is integrable on [2, +00). Define Slhp ^'^ ^ contour going from u to w = 2 in the lower half-plane 
bounded away from v = 1 and then from v = 2 to v = +00 along the real axis. Also define S^^p to be 
a contour going from v = +00 to w = 2 along the real axis, and then from w = 2 to = in the upper 
half-plane, bounded away from v — 1. By dominated convergence, we may pass to the limit y ^ — 1 in the 
integrand for l2'. 



/2(z/;Z,£,A^)=2-l/2+7/2e 



(83) 



y + 



g-*7r(l/2-7/2e)^-l/2+7/2£(^^ _ ^yE/e-f/^e-l 



+ / e-^^(l/2-7/2e)^-l/2+7/2.(^_^)»£/e-7/2e-l^^^^(^) 



T + 
'^UHP 



as 2/ — > —1. Deforming both contours so they lie on the negative real axis gives 



/2 = 2^/^+''/^" 



I sm TT 



(84) 



7 iE 

2£^~ 
2l/2+7/2e,^ 

r(i + f)r(i-i-f) 



(_„)-l/2+7/2e(_^ + l)«i^A-7/2s-l ^ ^(^) 

) 

(-«)-l/2+7/2e(_„ + ^y<B/e-7/2.-l + 



by the identity (|53|). Using the change of variables w = w/(u' — 1) and the identity (551 



(_„)-l/2+7/2e(_^^l).£/e-7/2.-l^„^ / ^-l/2+7/2e(i _ ^)-l/2-^E/e 



(85) 

Also, as y — > — 1, 
(86) 



r(l + 2^)rG-f) 
r(i + i-f) 



Using equations (76), (61 1, (63 1, and (86 1 and the factorial identity 
(87) 



(1 + y)'^l'^^ = 2'E|2e^^Ex|e^^ ^ ^^yy_ 

\ the factorial ide: 

r(i + z) = zr(z) 



therefore gives equation ( 74 ) . 



□ 



Proposition 2.6. The eigenvalues in the upper half of the complex z-plane are 

1. (Antikink) z = + 1 — ii)i, 

2. (Kink- antikink pairs) z ~ e^i with 6 G M satisfying cosh(&) = + 1 — 2ne for n G Z+ smc/i f/iaf 
1 < n< {^/JFTT- l)/2s, 

3. (Breathers) z — e*^ wii/i sin(^?) = fi^ + I — 2ne for n G Z+ smc/i i/iaf ( ^/i^ + 1 — l)/2e < n < 

+ l/2e. 



Proof. By general scattering theory. 

and 



5*22 (z) has an analytic extension from the real line into the upper half- 
by definition, the eigenvalues are the zeros of this analytic continuation in 



plane (see Theorem A. 7 1 

the open upper half-plane. For the special case of the initial data ([2]), the analyticity of 522(2:) can be seen 
from the explicit formula (74 1. Indeed, T{z) has no zeros and simple poles at 2; G — Z+. It follows that 
r(i — ^) and r(^ — ^) have no poles for z in the upper half-plane. Therefore 
exactly {\/ ij? + I — (case 1) and the poles of r(l — 57 — ^) (cases 2 and 3). 



the zeros of 6*22(2) are 
□ 



Proposition 2.7. The coefficient Si2{z) — £'12(2; £, /i) at t — is given by 

r(i-f)r(i + f) 



<S'l2(2) 



2 G 



Proof. Assume z G M\{0}. Using equations (75 1 and m2\ 



(89) 



Su{z) = flf.l''i] = lim det[jrJ+]=- lim (J+e^'^^^n 



X^ — OQ 



lim 2-^^/-'C2 / (s-l)-l/2+^/2^(5+l)-l/2+7/2e(^_yyiJ/e-7/2e^5. 



To analyze in the limit x — > —00, we begin by deforming S+ to the contour U S^"*", where E^"*" is 
the contour running in a vertical line from + i to —1 — i, and S^"*" is a horseshoe-shaped contour running 
from —1 — z to 2 to —1 -t- i, staying bounded away from s — 1, y, and —1. See figure 2(b) For s G 



(90) Is - 2/1-^/2- = \{S +l)-{y+ l)r^/2^ = {\s+l\^ + {y+ l)2)-7/4e < |, + ^-f/2e 

because 7/e > 0. Then 



(s-1) 2 + 2e(s-fl) 2 + 2e(s-y)e 2e 



< |s_l|-5+2t|5+l|-5 + A 



(91) 



(s-y)--'- 



for some constant > 0, and again, since 7/e > 0, the integrand is integrable on S^+. Also, the integrand 
is bounded and therefore integrable on E^^. Thus, by dominated convergence. 



(92) 



512(2) = 2-^^/^C2 / (5 - l)-l/2+^/2^(s + l)-l/2+^B/e 



Next, deform E+ to the contour running from —1 to 1 on the real axis along the lower edge of the branch 
cut for (s — l)^i/2+7/2£ then from 1 to —1 on the real axis along the upper edge of the branch cut. 

13 



Using the change of variables s — 2w — 1 and equation ( 55 1 , 



/E+ 



(s-l)-5 + K(s+l)-2+^ds: 



= -^^(-H2t)(l _ 5)-5+2t(l + + ds 



(l_s)-2 + i(l + s)- 



(93) 



:2isin TT 



2 2e 



(l_s)-5 + 4(l + s)-3 + ¥ds 



-1 



X + lE . . / / 1 7 

E ^ e I sm TT 

2 



.21+* + ^ 



r^2e^ e ism\ TT 



2e 
2 2£ 



(1 — W) 2 + 2e 



Therefore, 
(94) 



5l2(z) 



1 /I iE\ (\ iE 



rQ + ^)r(i + f) 
r(i + i + f) 



2 2e 



which completes the proof after the use of the reflection identity ( 53 1 



□ 



Proposition |2.7| gives immediately 



Proposition 2.8. Si2{z) = for e = £n{i-j) (see where N € Z+ 



The significance of this result, combined with Proposition A. 8 relating S2i{z) to 5*12(2;), is that e ~ 
7,7/3,7/5, . . . gives a sequence of values of e tending to zero for which the reflection coefficient p(z) :— 
S2i{z) / S22{z) is identically zero, and thus the scattering data are reflectionless and the corresponding so- 
lution of the Cauchy problem can be constructed from discrete spectral data only. The inverse-scattering 
transform may then be carried out more or less explicitly, a calculation we will perform in Section [3] 

Together, the formulae (74) and (88 1 show that p{z) = S2i{z) / S22{z) admits, in this special case 
of the initial data ([2]), a meromorphic continuation into the upper half-plane (generally S2i{z) admits no 
continuation of any kind from the real axis z G K). The meromorphic continuation of p{z) to the upper 
half-plane that is available in this case will have poles not only at the zeros of 5*22(2) (these are, by definition, 
the eigenvalues), but also at the poles of 5'2i(z). These latter poles are those of F (| + ^) (again using 
Proposition A. 8 to relate 5*21 to 5*12); in the upper half-plane these are: 

• (Imaginary axis) z = e^i with 6 € M satisfying cosh(6) = (2n-|-l)e for n E Z+ such that n > (1— e) /2e. 



(Unit circle) z 



with sin(0) = (2n + l)e for n G Z+ such that n < (1 — e)/2e. 



These "phantom poles" (poles of 521(2;) in the upper half-plane) are not (necessarily) eigenvalues. However, 
when the reflection coefficient is nonzero, they will affect deformations of the Riemann-Hilbert problem that 
are used in asymptotic analysis. See the discussion at the end of Appendix [A] 

Next we calculate the proportionality constants {rjn}, defined by ji {x; z) — 77„j2 {x; z) when z is an 
eigenvalue, for the reflectionless cases e — SNifJ-)- 



Proposition 2.9. Let e — eAr(/i) (see ^) where N G Z"*". Let z he an eigenvalue in the closed first quadrant 
and set n — ('j + 2iE(z))/2e (note n G N). Then the corresponding proportionality constant is ry„ = (— l)""""^. 



Proof. From equation (180), Aji {x; z) — 77nAj2 {x; z), where z is the above eigenvalue indexed by n, and rjn 
is its associated proportionality constant. Thus 



(95) 



j2 (2;; 2:) = Vnit{x;z). 



The second entry gives in particular J22{x;z) — irinJ22{x',z). We evaluate equations (62) and (73) at e 
£n{p)- Note that 



(96) 



T:=„_2iV-l and -^- + ^^N. 
£ 2e 2 2e 
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With these substitutions, equation (62 1 takes the form 

(97) J+(x;z,£^(M),/i) = C2(l-2/2)-("-^^-i)/2 / {s+iris-iris-yr-'^'-Us. 

Since iV is a nonnegative integer, we may deform the contour S"*" away from s — ~l and s = 1 to a small 
circle around s — y. Thus 

(98) J+{x■,z,e^{^,),^,) = C2{l-yT^-'''-'^/' f (s + l^s - ifis - yr-'^'-Us. 



Likewise, equation (73 1 becomes 

{x; Z, e^(/i), m) = C2(1 - y2)-(n-2Ar-l)/2 I _ ^)Ar(^ _ _ ^)„-2^-l 

JT.- 

(^^9) ^ (-1)«- ^2(1 - y2)-(„-2Ar-l)/2 /■ ^ ^)JV(^ „ ^)Ar(^ _ 2;)"-2A^-l ds 

= (-l)""V2+2(x;Z,£Ar(^),/i), 

and so rjn ~ (^1)"^^ by comparison with equation ( [98| . □ 

The results for the scattering data are summarized in Theorem |2.10[ 

Theorem 2.10. The scattering data for the sine-Gordon equation ([T|) at t — with initial condition ^ are 
as follows. 

and 

(102) Sn{z) = S22{-z), S2i{z) = -S,2{-z) ■ 

The eigenvalues in the upper half-plane are 

1. (Antikink) z = (-yZ/i^ + 1 — fi)i, 

2. (Kink- antikink pairs) z = e^'i on f/ie imaginary axis with 6 G M satisfying cosh(fe) — -yZ/i^ + 1 — 2n£ 
/or eac/i n G Z+ smc/i f/iai I < n < (-yZ/i^ + 1 — l)/2e, 

3. (Breathers) z — e*^ on f/ie unif circle with sin(0) = + 1 — 2ne for each n G Z+ satisfying 
(v/m^TI - l)/2e < n < Va*^ + l/2e. 

T/ie eigenvalues are generically (with respect to fi (li M. and e > 0) all simple. The scattering data are 
reflectionless for e = ejv(M) where N G Z+. In the reflectionless cases, the proportionality constants are 

e 



rjn = (—1)" ^, where n = (7 + 2iE{z))/2e G N for any eigenvalue z = z„ in the upper half-plane, and th 
modified proportionality constants are = rjn/ S22{zn) where 

{N\f 



n = 

(IIRI S' (r. l-i2i(7 + f)(2W)!' 

^47^ (2iV-n)! z„ + (7 + Ai)iV 4/ 

Remark. More generally, for reflectionless potentials, 

(104) S'22{Zn) = n - ^fc) n(^« - ^^)"'' 

where the product runs over all eigenvalues (presumed simple) in the upper half-plane. From the point of 
view of the numerical inverse-scattering method used in this paper, the specialized formulae ( 103[ ) (which are 
adapted to the initial data ([2])) are especially useful because many of the factors in (1041 involve differences 
of nearly-equal numbers which lead to loss of accuracy in finite precision arithmetic, while the products in 
(1041 have been converted into products of integers in (1031 that can be evaluated with exact arithmetic. 



On the other hand, the more general formula (1041 leads to a representation of the modified proportionality 
constants {c^} as residues of a meromorphic function, and such a representation is useful in the context of 
deformations introduced to study the asymptotic (N — s- oo) behavior of the meromorphic Riemann-Hilbert 
problem of reflectionless inverse scattering [12] . 

According to Theorem 2.10 the eigenvalues for the sine-Gordon problem with initial data ^ lie on 
the imaginary axis and the unit circle. Those on the positive imaginary axis come in pairs symmetric with 
respect to reflection through the unit circle, except for a single distinguished eigenvalue at (v'm^ + 1 — fJ-)i- 
This eigenvalue contributes the net topological charge —1 of the solution u. The eigenvalues on the unit 
circle have imaginary parts that are equally spaced. The plots in figure [3] illustrate the location of eigenvalues 




-2 2 -2 2 -2 2 

(a.) n = 1, N = 4: (h) fj, = 1, N = 8 (c) ^ = 1, iV = 16 



Figure 3. Eigenvalues for refiectionless cases s — Sn{^) (see ([9|) with fixed /i = 1 and 
varying A'^. The gray circle is \z\ = 1. 

for /i=lase>Ois varied. If /i is varied as a parameter, pairs of eigenvalues on the unit circle in the 
upper half-plane corresponding to a single breather may collide ai z = i and bifurcate off onto the imaginary 
axis, forming a kink-antikink pair (see figure El). Note that this bifurcation preserves the total topological 
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Figure 4. Bifurcation of a breather quartet into a kink-antikink pair in refiectionless cases 
with e = £Ar(/i) (see (l9|) holding N fixed and varying ji. The gray circle is \z\ — 1. 



charge of u. Note also that for e and /i such that (7 — l)/2e e {1, 2, 3, . . . }, there exist double eigenvalues 
ai z — ±i. The existence of eigenvalues with algebraic multiplicity greater than one is worth noting. For 
instance, the self-adjoint Schrodinger eigenvalue problem associated with the Korteweg-de Vries equation 
admits only simple eigenvalues. 
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3. Inverse-scattering for reflectionless initial data 



We now reconstruct the matrix G(z) — G(z; x, t, e, /i) (see equation (199l) corresponding to the 
specific initial conditions Q from the exact scattering data given in Theorem |2.10 in the reflectionless case 
when p{z) := S2i{z) / S22{z) = 0. We therefore fix iV G Z+ and set e — eAr(/i) (see (9])). We also assume the 
condition that (7 — l)/2e ^ {1, 2, 3, . . . } so all the eigenvalues (poles of G(z)) are simple. Define 



(105) 



M 



1-1) 



as the number of kink-antikink eigenvalue pairs. Label the eigenvalues in the closed first quadrant as follows: 
zo = {\/Jj^ 



1- 20 = (\/m^ + 1 - m)* 

2. Zn — ie^ and Zn = ie^^ for b = arccosh(y^/? 



A0 



for 



1 - 2n£), I <n< M, ne 
~ 2ne), M + 1 <n < N, n e Z+. 



Note that for the purely imaginary eigenvalues in case 2, the meaning of our notation is that |z„| > 1 while 
[znl < 1. In this case, — z„ and — are also eigenvalues. In case 3 on the unit circle, if z„ is an eigenvalue 
then — z*, —Zn, and z* are also eigenvalues. 

3.1. Numerical linear algebra algorithm for reflectionless potentials. We use the conditions of the 
Riemann-Hilbert problem of inverse scattering (see Appendix [a|) to determine the matrix G(z) for {x, t) S . 
In any reflectionless inverse-scattering problem, G(z) has no jump discontinuity across the real z-axis, and 
is therefore a meromorphic function with poles at the eigenvalues. With the assumption that the poles are 
simple, we may therefore expand G(z) in partial fractions as 
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G(z) = 
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The superscripts on the constant matrices K 



u 

n 1 



K^, K^, and associated with the kink or antikink 
eigenvalues indicate if the associated eigenvalue is in the upper or lower half-plane, and the superscripts on 
the constant matrices B^, B^^, B^^^, and B^^ associ ated with the breather eigenvalues indicate the quadrant 
of the associated eigenvalue. The residue conditions (2011 show immediately that the second columns of K: 

Write 
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n ) 



K^, B^j, and B,^ and the flrst columns of K^, K^, Bff , and B 
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vanish for all 
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for I < n < M and n — when applicable, and 
(108) B^„ = ' 

for M -|- 1 < 71 < A^. From the symmetries of equation (|18l| it follows that 
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(109) 



< 71 < M 

1 < 71 < M 



= -B,f * - -a2Bfa2 - CT2B7V2, M + 1 < 71 < iV. 



Note that the elements of K^, Kj^, K^, and are all imaginary. These symmetries show that the elements 
of the second row of G(z) can be expressed in terms of the elements of the first row, so to build G(z) it is 
sufficient to find the first row. Moreover, according to Proposition |A.16[ the potential u may be recovered 
from the first row of G(z), and in terms of the partial- fraction expansion (106 1 this results in the formulae 
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sin(u) 



(111) 




Recall that each eigenvalue z„ in the upper half-plane has an associated modified proportionality constant 
c„ which depends parametrically on x and t via an exponential factor. We denote by {c„} those modified 
proportionality constants associated with the eigenvalues labeled {z„}. Define the vectors 



(112) 
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■ • J cjv, ci, . . 
, . . , zjv, zi, . . 
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Applying the residue conditions ( |201 1 to the partial fraction expansion (106) yields a linear inhomogeneous 
system for p'-'^^ and p*^^^: 



(113) 



F 



F 

l2Ar+l 



p(l)- 




02Ar+l 


_p(2) 




a 



where Fi, 



0.1 



Here 02Ar+i is the vector of zeros of length 2N + 1 and l2Ar+i is the 2N + 1 by 27V + 1 identity matrix. The 
{x, t)-dependence of the coefficient matrix and the right-hand side of this linear system enters only through 
the modified proportionality constants making up the vector a. Eliminating p*^^'' using the first (block) row 
gives p'-"'^' = — Fp(^', and the resulting system for p'^^ is 

(114) (I + F2)p(2) =a. 

With the explicit use of the discrete scattering data, all of the entries of a and F are known functions 



of X and t. Thus, for any choice of x = a;o and t — ig, the system (113) can be solved numerically, giving 



(via equations (110) and (111 )) the value of u(a;o,io) independently of the value of u{x,t) at any other x or 
t values. 

3.2. Numerical results. Here we apply this procedure to study the solution of the Cauchy problem for 
the sine-Gordon equation ([T|) subject to the initial data (|2| for various values of the parameters fi and e that 
make the scattering data refiectionless (so that e = eAr(^) for some integer N). We are especially interested 
in the limit of large N, as this corresponds to the semiclassical limit of e J, 0. 

For large N, the system (114) is poorly conditioned, and it is therefore necessary to compute I + F^ 
and a with high precision at a given pair {x,t) G to find p'^-* and hence u{x,t) to even a few decimal 
places of accuracy. For instance, for N — 16, ^ — 1, x — 0, and t — 5, the condition number of I + F^ 
is approximately 3.5 x 10^^^, and it is necessary to use approximately 125-135 digit precision to accurately 
compute u. 

We first study the special case of /i = 0. Figure |5] shows plots of the square region —2.5 < x < 2.5 
and < t < 5 with different colors indicating different values of cos(u), with different plots corresponding 
to different values of N varying between N = and N = 16 (e = eAr(O) between 1 and 1/33). Lighter 
colors correspond to values of cos(u) closer to 1 and darker colors correspond to values of cos(u) closer 
to —1. The solutions of the sine-Gordon equation ([T]) illustrated in these plots consist of a "nonlinear 
superposition" of N breathers and one antikink. As each associated eigenvalue lies exactly on the unit circle 
in the z-plane, the velocity of each of these soliton components, when considered in absence of the others, is 
exactly zero. In this sense, the solution may be considered as a zero- velocity bound state of N breathers and 
one antikink. The most interesting phenomena are associated with the semiclassical limit £ | equivalent 
to letting N (the number of breathers) tend to infinity. In this limit, the plots suggest the asymptotic 
emergence of a fixed caustic curve t = t(x) in the space-time plane separating regions containing different 
kinds of oscillatory behavior. Indeed, for sufficiently large (that is, outside of the caustic), one observes 
roll patterns characteristic of single-phase traveling waves. The latter are simply the exact solutions of the 
sine-Gordon equation obtained by substituting into ([ij the traveling-wave ansatz u{x,t) — f{{kx ~ ujt)/e), 
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(e) fj. = 0, N = 8 {e = 1/17) ({) ^ = 0, N = 16 (e = 1/33) 

Figure 5. Plots of cos(w) for fixed fi — and varying N with —2.5 < a; < 2.5 (horizontal 
axis) and < t < 5 (vertical axis). 



resulting in the ordinary differential equation 

(115) ^u:^-k')f^+smif)^0.. 



^{kx — cut) . 



Here k is the wavenumber and uj is the frequency of the traveling wave, and the waves appearing as the roll 
patterns in figure [s] correspond to phase velocities c :— uj/k with |c| > 1 (which makes (1151 a time-scaled 
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version of the simple pendulum equation)]^ The periodic solutions of (115) are expressed in terms of elliptic 



functions, and therefore we say that the roll patterns in figure [5] correspond to modulated waves of genus 
G = 1. In the context of the phase portrait of the simple pendulum, the roll-pattern oscillations outside of 
the central region enclosed by the caustic curve correspond to librational motions of the pendula, i.e. orbits 
inside the separatrix. The sine-Gordon equation ([l]) also has families of exact solutions associated with 
hyperelliptic Riemann surfaces of arbitrarily large genus G, and these solutions are represented in the form 
u(x, t) — f(9i/e, . . . , Qg/s) where On — knX — ujnt and where / is a multiperiodic function of period 2t: in 
each of its G arguments. In the case G > 1, m is no longer a traveling wave, but rather is a multiphase wave. 
Reasoning by analogy with understood semiclassical limits of other integrable equations, we may expect that 
the more complicated oscillations evident in the plots of figure [s] for t > t(x) (that is, inside of the caustic 
curve) are modulated multiphase waves for some G > 1. Finally, we note that the caustic curve t = t(x) 
appears to originate from the point x = t = 0. As the velocity Ut is zero at t = and the pendulum angle is 
u — —IT at X ~ t ^ 0, the point x = is the unique point in the initial data corresponding to a point on the 
separatrix of the phase portrait of the simple pendulum. 

The evolution of the initial data ([2| for 7^ is depicted in the plots shown in figure |6] These plots 
are analogous to those in figure |5] except that we set — 1 and considered e = eAr(l). The main effect of 
nonzero fi on the discrete spectrum is to include, along with the quartets of eigenvalues that correspond to 
breathers, an asymptotically (in the limit N 00) nonzero fraction of eigenvalues on the imaginary axis 
that correspond to kinks and antikinks. The velocities of the kinks and antikinks asymptotically fill out 
the entire range of values (— |^|/-\/^I^~-i-T, ImI/V + !)• There is always one more antikink than there are 
kinks, and the "excess" antikink (corresponding to the eigenvalue Zq = — l)i) carries the topological 
charge. This excess antikink always moves to the right (this is a consequence of /i > 0, it turns out), and the 
kink-antikink pairs corresponding to the other eigenvalues on the imaginary axis are shed periodically in time 
and move to the left and right. As e J, 0, the outermost kinks or antikinks form a caustic curve separating the 
modulated single-phase waves outside from a region of the space-time containing the kink/antikink trains. 
As these trains propagate outwards over a field of modulated single-phase waves, it seems reasonable to 
suppose that the pattern in this part of the space-time would be described by a modulated multiphase wave 
of genus G — 2 that may be viewed as a nonlinear superposition of the single-phase waves (G — 1) and a 



kink or antikink train (also G = 1, although via orbits of (115) in the case w^/fc^ < 1 that lie outside of 
the separatrix). That the antikinks are moving to the right while the kinks are moving to the left (for these 
plots corresponding to /i = 1 > 0) can be seen from a plot of u itself reconstructed from its sine and cosine 
subject to the boundary condition u{—oo,t) = as shown in figure [Tj This plot corresponds to a horizontal 
slice of figure [6(f)] (or 8(c) below), and it is completely clear that the kinks occupy the left-hand portion of 



the plot (in which from figure [6(f)] we see that the waves are propagating to the left) while the antikinks 
occupy the right-hand portion (and by similar observations are propagating to the right). 

The caustic curve simultaneously emerges at t = from two (asymptotically) symmetric nonzero 
points X, and again these points admit an interpretation in terms of the separatrix of the simple pendulum 
equation (see below). Between these G = 2 regions there is a triangular region containing pure single- 
phase oscillations that persists for a time independent of N. In the context of the phase portrait of the 
simple pendulum, these oscillations correspond to rotational motions of the pendula, i.e. orbits outside the 
separatrix. The collision of the two G — 2 regions at the top of the triangular G = 1 (rotational) region 
results in a region containing more complicated oscillations that resembles the region inside the caustic 
curve for /i = as seen in figure [5] Note, however, that the oscillations occupying this central region may be 
expected to be even more complicated than those present for fj, = because there are many kinks/antikinks 
with very small velocities, and these will (if N is sufficiently large) begin to interfere with the bound state 
of breathers. Finally, note that, while for /i 7^ the exact solutions are not symmetric about x ^ 0, the 
asymptotic behavior evidently becomes symmetric as e | 0. 

The effect of varying /i can be seen from the plots shown in figure [8] Here, N is fixed at the value 
= 16 and /i is varied, with e = Sie{fi) holding to ensure a reflectionless potential. Note that the base 
of the triangular region of the space-time containing single-phase rotational oscillations appears to increase 
with /i. The ratio M/{N — M) of eigenvalue quartets corresponding to kink-antikink pairs to eigenvalue 



"^Of course, for these solutions of the hyperbohc sine-Gordon equation, the phase velocity exceeds the light speed of |c| = 1. 
In a sense, this fact does not contradict the hyperbolic nature of the equation, because the traveling wave is certainly not 
spatially localized, and moreover it has infinite energy. 
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(a) fj. = 1, N = {e = V2) 




(c) = 1, AT = 2 = V2/5) 




(b) /X = 1, Af = 1 = V2/3) 




(d) /X = 1, Af = 4 = V2/9) 




(e) /i = 1, TV = 8 (e = v^/It) 
Figure 6. Plots of cos(u) for fixed n — 1 and varying 

axis) and < t < 5 (vertical axis). 



(I) H = l, N = 16 (e = V2/33) 
N for -2.5 < x^< 2.5 (hdr 



(horizontal 



quartets corresponding to breathers (see ( 105[ )) also increases as /i increases for fixed N, an effect that is 
clearly visible in the plots of figure [8j 

The plots in figure [8] also contain annotation indicating our best guesses as to the values of x from 
which the primary caustic curve emerges ai t — 0. These x- values may be predicted by the following 
simple argument. Let us rewrite the sine-Gordon equation ([T]) as a perturbed simple pendulum equation in 
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Figure 7. Plot of u ior t ^ 2.5, -2.5 < a; < 2.5, iV = 16, /i 1 (e = ^2/33). 




(a) = 0, Af = 16 




-0.89 0.88 

(c) 11 = 1, N = 16 




-0.48 0.45 

(b) n = 1/2, AT = 16 




-1.47 1.44 

(d) 11 = 2, N = 16 



Figure 8. Plots of cos(u) for various n with —1.5 < a; < 1.5 (horizontal axis), < t < 5 
(vertical axis), and N = 16. 



first-order form: 
(116) 



du dv 

e— = V , e- 



dt " ' " dt '^Mu) + £^F{t;x) 

with forcing term F{t] x) :— Uxx{x, t). We think of u and v as the angle and angular velocity of a pendulum 
indexed by a parameter x. At the initial instant of time < = 0, the function u is smooth and independent of 
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e, so the perturbation term e^F{t;x) is very small, and one expects u{x,t) to evolve nearly independently 
for different values of x € M. This situation of independent pendulum motions might be expected to persist 
until u develops rapidly-varying features of characteristic length proportional to e, for in such a situation 
we would have d hence the perturbation term is no longer negligible compared with sin(u). 

Now, at any fixed time t, we may plot the phase points (u, v) in the phase portrait of the simple pendulum 
(that is, of (1161 with F = 0), and this data will appear as a curve parametrized by x. Figure [9] shows the 
initial data ([2| plotted parametrically in the phase portrait of the simple pendulum for /i = 0, /i = ±1/2, 
fj, = ±1, and fi = ±2 (blue curves). The separatrix for the simple pendulum equation is shown with red 




Figure 9. The initial data ^ plotted for /i 
phase portrait of the simple pendulum. 



0, ^ = ±1/2, fJ, — ±1, and n — ±2 in the 



curves. It is clear that each blue curve intersects the separatrix at exactly two points, and moreover, by 
unraveling the parametrization it is easy to see that these two points correspond to two distinct values of 
X. Near these values of x, there are pendula undergoing librational motions as well as pendula undergoing 
rotational motions. This is the scenario under which the most rapid amplification of the difference of angles 
u for neighboring pendula is to be expected. Therefore, we may make the prediction that the modulated 
single-phase ansatz should break down immediately at i = at exactly the two values of x at which the initial 
data meets the separatrix. These values of x are easily calculated. Indeed, the separatrix is given by the 
equation v — ±2cos(u/2), and the initial data satisfies cos(m/2) = tanh(x) and v = 2/isech(x). Therefore, 
the initial data curve (blue) intersects the separatrix (red) at values x for which 



(117) 



fi — ± sinh(a;) 



To confirm this reasoning, we took our best guesses for the x-values at which the phase transition occurs at 
t = as indicated on the plots in figure [8] and created a data set by combining these with the corresponding 
values of fi. The ordered pairs (x, fi) making up this data set are plotted with black dots in figure 10 along 
with the curves (1171 plotted in red. It is clear that this theory provides an accurate prediction of the points 



X from which the caustics emerge at time t — 0. While we have only given a comparison with the theory for 
initial conditions of the special form ^ , it seems reasonable that the principle should be the same for more 
general initial data. That is, one should locate the x-values at which the pair {u — f{x),v — g{x)) lies on 
the separatrix v — ±2cos(u/2) and expect complicated oscillations to emerge from these points for t > in 
the semiclassical limit. 

We have only computed solutions corresponding to the initial data (p| for /i > 0. That this is 
sufficient follows from a simple symmetry between /i and — /i. Indeed, write ([1| and ([2| in first-order form 



(118) 



eut = V 



evt — e^Uxx — sin(u) 
subject to the initial data 

(119) u{x, 0) = /(x) , v{x, 0) = 2n sech{x) , 
where sin(//2) = sech(a;), cos(//2) = tanh(a;), and consider the substitutions 

(120) V=-v, U = 2n-u, X = -x, T = t. 
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Figure 10. The relation (117) and the numerical data. The data was collected only for 
/i > 0, but by a natural symmetry (see below) we may include the point {—x, —^) whenever 
we measure the point {x,ii). This plot suggests that an important role is played by the 
separatrix of the simple pendulum in the development of caustics for the semiclassical sine- 
Gordon equation. 



Then the Cauchy problem for U and V consists of the first-order system 

eUt = V 

(121) n , , 

^ ' e^T = e^C/xx -sin(t/) 

subject to the initial data 

(122) U{X, 0) = F{X) , V{X, 0) = -2^ sech(X) , 

where sin(F/2) = sech(X), cos(F/2) — tanh(X). Therefore, U satisfies the sine-Gordon equation with initial 
data of the form ^ but with /i replaced with — /i. In terms of cos(u), replacing ji with —fi therefore simply 
amounts to replacing x with —a:. 

4. Concluding Remarks 

The main result of this paper is the exact calculation, via the theory of hypergeometric functions, of 
the scattering data for the noncharacteristic Cauchy problem for the semiclassical sine-Gordon equation ([ij 
subject to the initial data (|2|. That this calculation is valid for all sufficiently small e means that the formulae 
for the scattering data given in Theorem |2.10| may be used to formulate a corresponding inverse-scattering 
problem whose solution will give detailed information about the semiclassical limit of the sine-Gordon Cauchy 
problem. Moreover, since for each value of the parameter S M appearing in the initial data ([2| there exists 
a sequence {e — eM{lJ)}^^i of values of e tending to zero for which the scattering data are reflectionless, 
it is possible to approach the semiclassical limit in such a way that the inverse-scattering problem involves, 
for each N , only finite-dimensional linear algebra. As we have shown in Section [3] this fact makes it quite 
feasible to use numerical methods to solve the inverse-scattering problem for fairly large values of N and 
therefore study the semiclassical limit, at least in a qualitative sense. Our numerical reconstructions of 
the exact solutions of the Cauchy problem indeed reveal marvelous structures apparently emerging in the 
semiclassical limit. 

Needless to say, a study of the semiclassical limit based solely on numerics of the sort described in 
Section [3] has practical limitations. To study the semiclassical limit really requires allowing N to become 



arbitrarily large, and the system (1141 contains 2N + 1 equations and hence will ultimately become numer- 
ically intractable for sufficiently large N . This difficulty is compounded on the one hand by the fact that 
the condition numbers of the matrices involved grow rapidljj^ with N , and on the other by the necessity 
to use a grid spacing of order e to resolve the microstructure of the solution. In other words, to study 
the semiclassical limit in this way, an asymptotically badly-conditioned linear algebra problem in dimension 
proportional to N must be solved on a grid of approximately ~ iV^ values of {x, t) in a fixed-size region. 



"'One can see from the formula for the matrix elements of F ( |113| that F is proportional by the diagonal matrix diag(a) to a 
matrix of Cauchy/Hilbert type. The latter is the classic example given in textbooks on numerical analysis of an ill-conditioned 
matrix. 
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In our opinion, the main purpose of carrying out numerical experiments like those in Section [3] is 
to indicate phenomena that would be of interest to study rigorously by other (analytical) methods, and to 
motivate such a study. For example, figures [5] and [6] clearly indicate the existence of a limiting form of the 
0(1) scale macrostructure independent of e in the semiclassical limit. The (apparent) existence of caustic 
curves separating different types of oscillations requires a careful explanation, and such an explanation would 
be expected to also make asymptotically accurate predictions for the locations of the caustics. In integrable 
problems like the sine-Gordon equation, one expects the microstructure of oscillations in between the caustic 
curves to be described asymptotically by modulated exact multiphase solutions of the equation associated 
with Riemann surfaces of genus G. The modulation itself is expected to be described by slowly- varying (that 
is, independent of e) fields satisfying an appropriate system of quasilinear Whitham (modulation) equations. 
The sine-Gordon problem is quite different from other integrable problems for which the semiclassical limit 
has been investigated in that it has Whitham equations of both hyperbolic and elliptic type ^ |7j. To 
fully analyze these phenomena from the starting point of the scattering data we give in Theorem |2.10[ it 
is necessary to use very precise methods of asymptotic analysis for Riemann-Hilbert problems to find an 
asymptotic expansion for u{x,t]£, jj) valid as e | 0. Calculations of this sort, also in the discrete spectral 
setting (that is, reflectionless inverse-scattering as is available for this problem when e — ejv(M))i were carried 
out for the semiclassical focusing NLS equation (|7| for a general class of initial data in [T^]. 

True understanding of the semiclassical asymptotics of the Cauchy problem for the sine-Gordon 
equation ultimately requires generalizing the one-parameter family of initial data given by While the 
special initial data ^ is quite natural, satisfying the correct boundary conditions, and incorporating effects 
such as nontrivial topological charge and tunable (via the parameter jj) initial velocity, one may certainly 
pose the Cauchy problem for more general initial data /(•) and g{-) and ask for the corresponding asymptotic 
behavior of u{x, t;e) as e J, 0. One might hope that other initial conditions that are somehow close to ^ 
might correspond to scattering data and dynamical behavior of u(x,t;s) whose semiclassical asymptotics 
are similar to those of the exactly solvable case. This would indicate a kind of stability of the semiclassical 
limit. Furthermore, one may be interested in the semiclassical asymptotics corresponding to initial data that 
differ significantly from the special data ([2]), for example by having a topological charge that is different. 
Clearly, to begin to study more general initial data, it is necessary to find quantitative approximations of the 
corresponding scattering data. A first step towards this goal is to seek conditions on general initial data that 
force the eigenvalues to lie exactly on certain contours in the complex z-plane. For initial data satisfying such 
conditions, WKB analysis can be used to find a leading-order estimate in e for the scattering data, and with 
more work, the error of the estimate can be analyzed. For example, for the nonselfadjoint Zakharov-Shabat 
eigenvalue problem relevant to the focusing NLS equation ^ , Klaus and Shaw J4j showed that if the initial 
condition q(x, 0) is real and monomodal, then the discrete spectrum may only lie exactly on the imaginary 
axis. In the semiclassical setting, this is an exact result that holds for all e > 0. WKB calculations based on 
the Klaus-Shaw result were used in [12] to analyze certain so-called semiclassical soliton ensembles. As for 
sine-Gordon, Bronski and Johnson [2] recently found a result analogous to that of Klaus and Shaw, showing 
in particulaij^ that if g = and / is monotone with topological charge ±1, then the eigenvalues must lie 
exactly on the unit circle in the z-plane. A quantitative approach to the discrete spectrum for initial data of 
the Bronski- Johnson type would be to map it onto a perturbation of the specific initial data ^ with /i = 
(which is, of course, a special case of a Bronski-Johnson potential) using a Langer transformation, and to 
control the error introduced by the perturbation for small e. This will be carried out in future work. 

We conclude by drawing some comparisons between our results and those of Tovbis and Venakides 
PT] for the nonselfadjoint Zakharov-Shabat problem associated with the focusing NLS equation. The class 
of Tovbis- Venakides potentials (see Q) involves a parameter /i € M (/i = is the special case studied earlier 
by Satsuma and Yajima |19j). In fact, we chose to use the symbol /i for the parameter in ^ precisely 
because this parameter plays a similar role. One immediate observation is that only in the case /i = is the 
Tovbis- Venakides initial data of Klaus-Shaw type, and similarly only in the case = is the initial data (pi 



Actually, they showed more: if g = and sin(//2) is a Klaus-Shaw potential, then the discrete spectrum lies on the unit 
circle. If the maximum value of this potential is unity, then the topological charge is ±1 and / is monotone, but if the maximum 
value is smaller the topological charge is zero and / is monomodal. 
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of Bronski- Johnson type. Thus only for fi — is one guaranteed by general argumentsj^ that the discrete 
spectrum is confined to a special curve in the complex plane. Another observation is that the parameter 
fj, has a physical interpretation of velocity in both the Tovbis-Venakides family of potentials (because in 
the hydrodynamic variables for Schrodinger equations introduced long ago by Madelung, the velocity of 
the quantum-corrected fluid motion is expressed in terms of = A{x,t)e^^^^'*^^'^ by Sx{x,t), and for the 
Tovbis-Venakides potentials S is proportional to ^ at i = 0) and also in the family ^ of initial data for 
sine-Gordon (because the initial data is a solution of the advection equation with velocity /i/e as pointed out 
in the Introduction). However, as one important distinction, we note that the Tovbis-Venakides potentials 
have the possibility of being reflectionless for certain e only for ii — 0, while this possibility exists for the 
initial data Q for every ^ e M. 
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Appendix A. The Riemann-Hilbert Approach to Inverse Scattering for Sine-Gordon 

Our aim in this appendix is to present a completely self-contained theory of inverse-scattering for 
the sine-Gordon equation in laboratory coordinates. In particular, we show how to represent the solution of 
the Cauchy problem for the sine-Gordon equation with L^-Sobolev initial data (specifically, 1 — cos(u(a;, 0)), 
sin(u(a;, 0)), Ux{x,0), Uxxix,0), ut{x,0), Utx{x,0) € L^) in terms of the solution of a certain matrix-valued 
Riemann-Hilbert problem. To ensure that various quantities used in the inverse-scattering method are well 
defined with desirable properties for alH > 0, we rely on a theory of the well-posedness of the Cauchy problem 
that may be developed independently of any inverse scattering methodology. An outline of the relevant well- 
posedness theory is given in Appendix [b] in which we show that the class of i^-Sobolev potentials (in the 
sense defined above) is preserved for all i > under the evolution of the sine-Gordon equation. Many of 
the results to be described below have appeared in the literature in one form or another. For instance, the 
characterization of the Jost solutions assuming that 1 — cos(M(x,t)), sin(ii(a;, t)), Ux{x,t), Ut{x,t) g at 
time t appeared in Kaup |13j . and aspects of the Riemann-Hilbert approach to inverse scattering were worked 
out for initial data / and g in the Schwartz space by Zhou (26j and Cheng et al. [HIS]. The well-posedness 
theory we present in Appendix [B] appears to be a new contribution to the subject. 

The starting point for our analysis is the observation [13^ that the sine-Gordon equation (jlj is the 
compatibility condition for the Lax pair 



(123) 



4iewj, = Lw 



4£; + z-i(l -cos(u)) 
-z^^ sin(ii) + ie{ux + Ut) —AE — z^^{\ — cos(u)) 



— z sm{u) — i£{ux + Ut) 



(124) 



4zeWf = Bw 



AD 



^(1 — cos(w)) 
ieiux 



z ^ sin(M) — i£{ux + Ut) 
-4L> + z^i(l-cos(u)) 



z ^ sin(u) -I- i£{ux + Ut) 

with D{z) and E{z) given in (151. In other words, there exists a basis (determined, say, by specification of 



two linearly independent vectors w at a; = i = 0) of simultaneous solutions of (1231 and (124) if and only if 
u = u{x, t) is a solution of the sine-Gordon equation ([TJ. 



The Lax pair (124l"(123l appears to have a singularity at z = 0. However, it is possible to use a 
gauge transformation to move the singularity from z = to z — oo and in this way analysis for large z can 
be continued to appropriate sets with limit point z = 0. This gauge transformation will play an important 
role in our analysis. 



^It turns out that for /i 7^ the eigenvalues of the Tovbis-Venakides potentials (when they exist) lie exactly on the imaginary 
axis nonetheless. However, for the initial data ^ any nonzero value of /i immediately introduces eigenvalues that are not confined 
to the unit circle. 
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A.l. Jost solutions of the scattering problem. We now attempt to define the Jos t sol utions J (x) = 
as the fundamental solution matrices of the eigenvalue equation ( 123 1 normalized as 

-iEx/e 



3^{x; z, t, e) for z e 



(125) 



3 (x) 





'^ — lEx/s 









We denote the columns of J (x) as 
(126) 



0(1) 
0(1) 

J"(x) =: [j^(a;),j^(2:)] 



as X 



as X —> —oo. 



The issue at hand is to determine whether these conditions uniquely determine J (x) when z is a real 
number, and then to further determine what can be said for complex z. 



To begin, we rewrite ( 123 1 in the form 



(127) 
with 



(128) 



Q{x;z,e) := 



z-\l - cos(/)) -z-i sin(/) - i{ef + g) 

^z-^sm{f)+i{ef + g) -z-i(l - cos(/)) / 

Here / is the value of u, and g is that of eut at some fixed time t. The purpose of this decomposition is to 
separate the part of the coefficient matrix that decays (in a certain sense) as x —^ ±oo (Q) from a constant 
term (AEas). Defining matrices 

(129) M'^= [mf,m^] J'^e*■^^'"^/^ 

or equivalently in terms of the columns. 



(130) 



m; 



Ji e 



lEx/e 



mi = Ji e 



lEx/e 



rri — i p^'^^^/'^ 



one may easily translate the differential equation ( 127 1 and boundary conditions ( 125 1 for z e E into integral 
equations for the matrices M (x; z, e): 

(131) M+(x) = 

(132) MTix) = 



1 


0' 


1 r+°^ 





1 




1 


0' 


1 r 





1 





~iE{x-y)a-3/e 



Q{y)M {y)e 



dy, 



-iE{x~y)a-3/e 



Q(y)M (y)e 



iE{x-y)a3/£ 



dy. 



While these integral equations are formulated to correspond to (1231 and (1251 for z g M, we may also 
consider them for complex z. Proposition A.l shows that the columns fn^{x) and ni^(x) are well-defined 



by (131) and (1321 respectively as long as 3(z) < 0, and moreover for each a:: e M they are analytic for 
3(z) < 0, and continuous in the closed lower half z-plane for z bounded away from z — 0. Then Proposition 
|A.2| uses an alternate gauge to extend continuity to small z. 



Proposition A.l. Suppose 1 — cos(/), sin(/), ef' + g E L^. // 3(z) < 0, then the first column of (1311 
and the second column of (132) uniquely define solutions ni'^{x; z) and in^(a;; z). These functions are, for 
each x £ K, analytic for 3(z) < and continuous for z g {|z| > (5} H {3(z) < 0} for each (5 > 0. 



Proof. The function ni'^{x;z) is constructed from equation (131) via an iterative argument. Define the O"^ 

e the n**^ iterate inductiv 

/ Ki (y;a:,z)m+„_i(2/)dy 

J X 



iterate for as Yn^^{x) \— [1,0]'''. Then define the n*^ iterate inductively by 



(133) 
with 

(134) K+(2/)=K+(y;x,z) 



m7„(a:) 



4ie 



z-i(l -cos(/)) -z-isin(/)-i(e/' + g) 

(-Z-1 sin(/) + i{ef + g)) e^^C^-y)-^/^ -z-i(l - cos(/))e2*(^-2')^/^ 
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Here /, and g are functions of y. It follows that 



(135) 



+ (-!)' 



-\-OC /' + 00 



+00 



K^(a;i) • • •K^(a;„) 



dxn ■ ■ ■ dxi. 



If the sequence {m^„(a;)} converges, then m^(x) will be defined as its limit, which clearly has the form of 
an infinite series. 

Consider the n**^ term in this series. Let |jv|| :— \vi\ + \v2\ be the £^ vector norm and ||M|| — 
max(|Afii| + |M2i|,|Mi2| + IM22I) be the induced matrix norm. The key observation is that (because 
Q{E) < for 9(z) < 0) the assumption 3(z) < implies that if y > x then ||K^(y)|| is bounded by a hnear 
combination of 1 — cos(/), | sin(/)|, and \ef' + g\ with constant coefficients independent of y and uniformly 
bounded for |z| > (5 > 0. Therefore whenever 5(2) < with z ^ we may define a function in L°°(M.) by 

f +00 



(136) 



\Kiy)\\dy. 



Furthermore, is uniformly bounded in z for z G {|z| > ^5} H {^(-s:) < 0} for every ^ > 0. Then 



+00 r-\-oo 



+00 



< 



(137) 



+ 00 /' + 00 



1 



+00 



dxn • ■ ■ dxi 



l|Ki(xi)||-- 

dv{Xn) ■ 



(a;„)|| dXn 
■ dv{xi) 



■ dxi 



iy(a;)" 



It follows that the partial sums m^„(x) are majorized by those of an exponential series, and so the sequence of 
partial sums converges and the limit furnishes the unique solution of the first column of the integral equation 
( 131 1. By uniformity of the convergence, analyticity for 9(z) < and continuity for z € {|z| > 5}n{9(z) < 0} 
for each S > extend from the partial sums to the limit in'l{x; z). We also have the estimate 



(138) 



||m+(^)|| < e--^^) < exp (^^^ W^tivWv 



< oo. 



which is uniform for 5 > 0. The argument for m.2 {x; z) is similar. 



□ 



The argument in Proposition A.l fails for z near z = because of the coefficient z ^ in the matrix entries 
of K]^ (y). The use of an alternate gauge, which we call the zero gauge, circumvents this problem. We define 



a new set of functions in terms of the Jost solutions J {x) by 



(139) 



J±(x) 



cos 



2 



sm 



2 



j*(x), 



with columns 3^{x) =: [j^(a;), j^(a:)]. Note that this gauge transformation can be interpreted as a rotation 
of the Jost solution column vectors by an angle —f{x)/2. It follows by direct calculation that the gauge- 
transformed matrices [x] satisfy the modified eigenvalue equation 

(140) 4i£j^ = (4^;ct3 + Q) 

where 

-z(l - cos(/)) -z sin(/) + i{ef' - g) 
-zsinif) ~ i{ef' - g) z(l-cos(/)) 
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(141) 



Q 



lim u{x, 0) 

X — *+oo 



2nn 



Assuming the boundary conditions 
(142) lim u{x,0)^0 and 

X — > — oo 

hold in a suitable sense, the required behavior of J_^{x) as x ^ ±00 is derived from (1391 and (1251: 





(143) 



J+(x) 
J-(x) 



(^_l'jn+l^-iEx/e 







-iEx/e 







^iEx/e 



(__l)"+lg*-Ex/s 

+ 0(1) as a; — 



+ 0(1) as x ^ +00 for z G 
—00 for z e M. 



Analogous to equation (1291, define 
(144) 



M^{x) := J±(a;)e*-^^'^^/= 



with columns M (x) =: [mf (x), xnf(x)]. It follows by integrating (1401 using the boundary conditions 
( 143 1 that (x) satisfy the integral equations 

f +00 



(-1)" 




(145) M^(a;) = 

(146) M-{x) 



1 











(- 


1)"+1 




"1 


0' 


1 







1 


+ — 

4te 


I 



1 



-iE{x-y)a3/e 



Q(y)M+(y)' 



JE{x~y)a3/e 



dy, 



-iE(x-y)(73/e 



Q(y)M-(y)e 



,iE{x-y)a3/s 



dy. 



Now these modified integral equations for the gauge-transformed solutions M (x) are used to show that the 
columns of M (x) are continuous in a neighborhood of z = in appropriate half-planes. 

Proposition A. 2. Suppose 1 — cos(/), sin(/), ef — g E . Then for each x E M. and for each S > 0, 
m]^(x; z) and m2"(x; z) are continuous functions of z in the region z G {\z\ < 6} O {^(2) < 0}. 

Proof. The function m^(x ) is constructed iteratively from equation (|l45[), similar to the construction of 



nL[(x) in Proposition A.l 

(147) 
where 

(148) K+(y;x,z) 



Set m+o(x) := [(-l)"+\ 0]"^. Define the 



iterate by 



m; 













Ki{y]x,z)m+^{y)dy , 



1 

4je 



-z(l-cos(/)) -zsiiiif) + i{ef' -g) 

(-zsin(/) - i{ef - g)) e^'^^-y')^/^ z(l - cos(/))e2^(^-2')^/^ 



Aside from the exponential factors e^**^^ y)E/e ^ everywhere that a factor of z ^ occurred in K^(y) there is 
in K^(y) a factor of z. This allows parallel analysis as in the proof of Proposition 



A.l 



to go through with 

the condition |z| > 5 replaced by the condition |z| < 5. Thus, the iterates converge and mj"(x: z) is analytic 
in the lower half z-plane and continuous in the closed lower h alf z- plane for bounded z. A similar argument 
works for (x; z) as well. Finally, the gauge transformation ( 139 1 is independ ent o f z and so does not affect 
the continuity, and therefore m^(x; z) and m^(x; z) as defined by Proposition 
the whole closed lower half-plane. 



A.l 



are in fact continuous in 
□ 



Together, Propositions 



A.l 



and 



A.2 



show m||*"(x; z) and rn^ix; z) are analytic in the lower half z- 
plane and continuous in the closed lower half z-plane. An analogous result holds for mj^(x; z) and m^(x; z) 
in the upper half-plane. 



Proposition A. 3. Suppose 1 — cos(/), sin(/), g E . //3(z) > 0, then the first column of (1321 and 
the second column of ( |131 1 uniquely define solutions rnj~(x; z) and m^(x; z). These functions are, for each 
X G M, analytic for 3(z) > and continuous for 3(z) > 0. 

Proof. The proof follows that of Propositions |A.1| a nd A.2 taking into account the sign of the real part of 
the exponential factors in equations (131) and (132). □ 



Propositions A.2 and A. 3 can be summarized as follows. 
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Theorem A. 4 (Kaup, [13 ). Suppose 1 — cos(/), sin(/), /', 5 e . Then rn.^ {x; z) andm2{x;z) are well- 
defined and for each x € M. are continuous for z e M and extend continuously to analytic functions in the 
upper half z -plane. Similarly, m.f{x] z) and fn^ix; z) are well-defined and for each a; e M are continuous for 
z GM. and extend continuously to analytic functions in the lower half z-plane. 

Next we establish a lemma showing that under the assumption of a little more smoothness of the 
potentials, the x derivatives of the columns of M (x) are, for each fixed a; G M, uniformly bounded in the 
appropriate closed half-planes for z. 

Lemma A. 5. Suppose 1 — cos(/), sin(/), /',/", g, 5' e . Then nll^{x; z) and m^^(a;; are uniformly 
bounded in x for each fixed z with 3(2;) < 0, having L°° norms that are uniformly bounded for all such z. 
Similarly, in^^{x; z) and m^^(a;; z) are uniformly bounded in x for each fixed z with > 0, having L°° 
norms that are uniformly bounded for all such z. 

Proof. We show the result for m^^. The proofs of the results for m^^, m^^, and mj^ are similar. Choose z 
such that < 0. From (131 1, satisfies the integral equation 



(149) 



xal{x; z) 



+00 



Ki {y;x,z)xa^{y;z) dy 



with K]^ defined by (134). We write the entries of Q and mf as 
(150) Q 



"Qii 


Q12 


, m+ =: 




Q21 


Q22. 


m+_ 



The first entry of ( 149 1 is 

(151) m+(x;z) = l- 
and differentiation in x gives 

(152) rn+,^{x;z) 



Aie 



+00 



Qu(2/; z)m+^{y; z) + Ol^iv, z)m+^{y; z) dy, 



^[Qn(a;; z)ml^{x; z) + (5^2(2;; z)m^^{x; z)]. 



Thus m^j^^(a;; z) is uniformly bounded in x with L°° norm uniformly bounded for z G {|z| > 5} n {S(z) < 0} 
because xa.^ and Q are (we see that f',g G L°° by noting f",g' G by assumption and applying the 
fundamental theorem of calculus) . The second entry of ( 149 1 is 



(153) 



■mti{.x;z) 



Aie 



+00 



g2.i=;(.-,)/s[Q^^(y. ^)jntM z) + Q22(y; z)m+ (y; z)] dy. 



Taking an x-derivative gives 



1 



^tixi^'^z) = —[Q2^{x;z)m+^{x;z) + Q22{x;z)m+^{x;z)] 



(154) 



1 

Aie 



+00 



2iE 



^''''^"""^/lQ2i(y; ^)m+i(2/; ^) + ^22(2/; ^)m+ (y; z)] dy 



4ie 



[Q2i{x\ z)m'l-^{x; z) 



z)Tr4^{x; z)] 



1 

Aie 



+00 



dy 



(e'^^^-yy-) [Q^^iy- z)mt,{y; z) + Q^-^iy- z)m+ (y; z)\ dy. 



Now (for y > X and in the indicated region of the z-plane) we have exp{2iE{x — y)/e),fn^-j^,fn2i G L°° . Also, 
Q ^ as a; — > ±00. To see this, note that the limiting value of f^ exists as a; — > ±00 because f^x G L^, and 
moreover since f^ G both limits must be zero. The same reasoning holds for 1 — cos(/), sin(/), and g. 
Therefore, 



(155) 



lim 1 — cos(/) — 0, lim sin(/) = 0, lim f^ = 0, lim 5 = 0. 

X — >±oo X — ^±00 X — ^±00 X — *±oo 
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Thus exp(2i£'(2; — y)/e)[(52i(y)^'^ii(y) + Q22(2/)'^2i(?/)] ^ as — > +00, and so integrating by parts and 
distributing the y-derivative gives 

(156) 

m+^(a;;z) 



1 

He 
1 

Aie 



+00 



+00 



lix) 



^2iE(x y)A[Q2iy(y; z)m+ (y; z) + Q^^{y; z)m+^y{y; z) + Q^^yiy; z)to+ (y; z)]dy 
J{y;x,z)m^^y{y;z)dy. 



Note that /(a;) is uniformly bounded in x with norm uniformly bounded for z G {|z| < 5} n {S(z) < 0} 
and that J ^ with norm uniformly bounded for z e {|z| < (5} n {^(z) < 0}. Therefore, by an iteration 
argument as in the proof of Proposition |A.1[ 

^+00 



(157) 



\Tr4ix\\L^ < II^IU- exp (^j \J{y)\dy^ < 00 



where the bound is uniform for z £ {\z\ < (5}n{3(z) < 0}. The uniform bound for z e {|z| > (5}n{5(2;) < 0} 
is shown similarly using the zero gauge defined by (1391. □ 

The assumption of additional smoothness of the potentials as above also provides limiting values of 
the columns of M {x; z) in various situations. 

Proposition A. 6. Suppose 1 — cos(/), sin(/), /',/", e L^. Then the columns o/M^(x;z) have the 

following limits in x and z: 

(158) 

lim ni^{x; z) = 

X — * + oo 
9(z)<0 

(159) 

lim mf{x \ z) - 

z — *oo 
Q(2)<0 



lim m.2 (x; z) 

X — * — 00 
O(z)<0 



lim m, (x; z) 

z — >oo 
3(^)<0 



lim (x; z) 

X — *— CO 

9(z)>0 



lim m, (x; z) 

z — >oo 
9(z)>0 



lim mj(a;;z) 



X — ^+00 

a(2)>o 



lim xa.t{x]z) 

z — >oo 

a(z)>o 



lim vol {x; z) — [—1) 

2— ►0 
9(2)<0 



n+l 



(160) 



2 

M 
2 



lim m2 [x; z) = 

Q(2)<0 



COS 



lim (a;; z) 

Z-+0 
9(z)>0 



sm 



2 
2 



lini m+{x;z) ^ (-1) 

z— >0 
9(z)>0 



n+l 



COS 



M 

2 

2 

2 

2 



Proof. We will prove the statements concerning {x; z) \ the proofs of the corresponding limits for {x; z), 
m^(x; z), and m^(x; z) are similar. 

We first establish the limit in x. Fix z E {\z\ > 6} {S(z) < 0} for some fixed (5 > 0. Consider 
the integral equation (1491 for ^^^{x; z). The product ni^ g as a function of y since \z\ > 6, since 
cos(/) — l,sin(/),/',(; € L^, and since rHj^ G for z € {|z| > (5}n{3(z) < 0}. Furthermore, K^m]^X[2;,+oo) 
tends to zero pointwise in y as a; — > +00. Therefore the limit for mf{x; z) as x — > +00 holds by dominated 
convergence for z g {|z| > 5} n {3(z) < 0}. The result for z G {|z| < 5} fl {3(z) < 0} holds by the same 
reasoning applied to the integral equation 

f +00 



(161) 



ni[{x; z) 



(_!)«+! 





/ K^(y;a;,z)m+(y;z)dy, 
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written in the zero gauge with given by ( 148 1, and the use of the gauge transformation ( 139 ) to go back 
to nif{x; z). 



Next consider the hmit of m.^{x;z) as z ^ oo for ^(z) < 0. The second entry of (149) may be 
written as 



(162) 



TO+ (a;;z) 



1 

iie 



+00 



dy \2iE 



^^E(.-v)/e\ [g^^(y. ,)jnUy. ^) + Q^^(y. (y; z)] dy. 



Integration by parts gives 



(163) 



4i(a;; z) ^^{Q^^(x\ z)mj^(a;; z) + Q2i{^\ z)mt^{x\ z)\ 



+00 1 A 

_,2.i=;(.-,)/s_[Q^^(y. ^)^+ (y. ^) + Q^^(y. z)mUy- z)\ dy. 



8E 



Since f',g E L°" (because f'^g' £ L^), we h ave Q ^^ , Q^^ E L°" for \z\ > 5. Since also and are 
uniformly bounded for ^(z) < by Theorem A. 4 the boundary term [Q2i™ii + Q22^2i\l^^ vanishes as 
z — > oo for 9(z) < (and hence as i? ^ oo). As for the integral term, Q211 Q2iyi Q22y G foi' 
\z\ > 5 and e^*^^^~^-'/^,m]^j^,m]^j^j^,m^j^,m^j^,y are uniformly bounded for y > x and 9(z) < 0. Therefore 
g2»iS(^-a)/e^/jy[Q_^^^+ _^ Q22wJi] G for z G {|z| > (5} n {5(z) < 0}. Since £; ^ cx) as z ^ cx), the 
integrand tends to zero pointwise in y almost everywhere as z — > 00 for 3(z) < 0. By dominated convergence, 

(164) lim m+ (a;;z) = 0. 

z — *oo 
9'(z)<0 



To analyze mj^(x;z) in the same limit, consider the integral equation (1511. The integrand is in for 
z e {|z| > (5} n {9(z) < 0} since Qn,Q^2 ^ for |z| > 5 and mti.mti E L°° for 3(z) < 0. In addition. 



the integrand tends to zero pointwise as z ^ 00 for 3(z) < since Quiy) tends to zero as z ^ 00 and ( 1641 
holds, while and f7iu,Qi2 G for 3(z) < 0. Thus, by dominated convergence, 



(165) 



(166) 



lim m{i{x\ z) = 1. 

Z— »C30 

9(z)<0 



Finally, we consider the asymptotic behavior in the limit z ^ 0. The statement that 



lim 1111(2::: z) = 

3(z)<0 



(-1)"+!' 




holds may be shown as above using the zero gauge. Then the limit of vclI{x\ z) as z — > for 3(z) < follows 
by inverting the gauge transformation with the help of (1391. □ 



Note that, from the asymptotic behavior of the columns of M {x\ z) in the limits x ±00 and the 



fact that (Abel's theorem) Wronskians of solutions of (1231 are independent of x, we have det(J {x; z)) = 1 

for a; e M and z G M. 



A. 2. Scattering data. The Jost solution matrices J (x; z) and J (a:; z) are both fundamental solution 



matrices of the same system ( 123 1 , so consequently the columns of J (a:; z) are necessarily linear combinations 
(with coefficients independent of a:) of those of J (a:;z). Therefore, there exists a matrix S(z) = S(z;t, e) 



such that 



(167) 



J (x; z) = J (a;;z)S(z), S(z;t,e) 



Sii{z;t,e) Si2{z;t,e) 
S2i{z;t,e) S22{z;t,e) 



z E 



The matrix S(z) is called the scattering matrix. The i-dependence of its elements comes from considering 
/ and g to depend parametrically on t (for example, if / = u and eg — ut come from a solution of the 
sine-Gordon equation ([ij). We will calculate this time dependence shortly (and in fact it will turn out that 
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the diagonal elements are independent of t). Using the fact that det(J {x; z)) — 1, we easily obtain the 
Wronskian formulae 



(168) 



S'ii(z) = det[ji (x; z), {x; z)] , 
S2i{z) = det[ji {x; z),j'^{x; z)] , 



Si2{z) = det[j2 (x; z),j^ {x; z)] 
522(2:) = det[ji {x; z), j^(a;; z)] 



These formulae, in conjunction with Theorem |A.4| and Proposition |A.6| provide a proof of the following. 

Lemma A. 7 (Kaup, [13]). Suppose 1 —cos(/), sin(/), /',/", 5' e . Then 822(2) is continuous for z e M 
and has a continuous extension into the upper half z -plane as an analytic function, while Sii{z) is continuous 
for z € M and has a continuous extension into the lower half z -plane as an analytic function. Moreover, 

(169) 



lim 5*11(2) = lim 6*22(2:) — 1 , 

z — *oo z — *oo 

9(z)<0 3'(z)>0 



and similarly 

(170) lim 5ii(z)= lim ^22(2) - (-1)"+' • 

z^Q z^Q 
a(2)<0 9'(2)>0 

Next we record several important symmetries of the scattering matrix. 

Proposition A. 8 (Kaup, '13]). For z G M, the elements of the scattering matrix are related by S*ii(z) = 
S22{~z) = 5*22(2)* and 512(2) = -52i(-2) = -521(2)*. 

Proof. Here it is essential that z e M so that both columns of J* (a;; 2) are simultaneously defined. In the 



eigenvalue equation (123), the coefficient matrix has the symmetry L(x; 2) = iT2L(a;; —z)a2. Therefore, 
(171) 



a23x{^'j ^z) = L(a;; 2)0-2 J"*" (2:; —2), 



and so J {x; z) — cr2J {x; ~z)C^ for some constant matrices C*. Write 
(172) M"^ {x;z)= 3^ [x; z)e'^'"''/^ 



Taking the limit as x — > ±00 and using Proposition A. 6 shows that C = 02 



Next, substituting the identity 



(173) 



3'^{x] z) = (72 J"*" (a;; —2)0-2 



8(2) = o-2S(-z)ct2, 



into equation (167) gives 
(174) 

which shows 5ii(2) — 822(^2) and 512(2) = — 52i(— 2). The matrix L also has the symmetry li(x; z) 
— L(a;, z*y that holds for all z G C, in particular 2 G M. Restricting to 2 G M, this implies 

3^{x; z)^ = —J^{x; z)^L{x; z). 



(175) 
Furthermore, 

(176) 



dx 



3 {x; zY 



-J^(x;z) ^ —3^ (x; z) ■ 3^ (x; z) ^ 

ax 



3 {x; z) ^h^x'jz), 

and so by comparison with (175) 3^{x; z)^^ = D*J^(a;; z)^^ for some constant matrices D*. Now 



(177) M'^{x;zY = e-'^'"''^'3'^{x;zY - e-*^^'^^/^D±j'^(a;; z)-i = e-'^="''^'-D'^e'^'"'^/'M'^ {x; z)-\ 



Again taking the limit as x — s- ±00 and using Proposition |A.6 | shows that 
(178) 



Substituting the identity 



3^ (x; zY = 3^ {x; z) ^ 



into equation (167) yields 
(179) 



S(z)t = S(2)-i, 



from which it follows that 5ii(2) = 522(2)* and 512(2) = —521(2)* 
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□ 



A.6 



By definition, the e igenvalues for the scattering problem ( 123 1 are the complex numbers z for which 
there is a solution of (1231 in 
to Proposition 



The Jost solution \^ {x\z) is defined for > and, according 



and the relation ( |130 1 between jj^ (a;; z) and nij^ (x\ z), (x\ z) decays exponentially to 

zero as a; — > — oo if and only if > 0. Similarly, the Jost solution (x; z) is defined for ^(z) > and 
decays exponentially to zero as a; — > +cx3 if and only if > 0. All other solutions blow up exponentially 
in these limits. Therefore, the eigenvalues z in the open upper half-plane are exactly those values of z for 



which j]^ (x\ z) is proportional to '^2 Recalling the representation (168) of 822(2) as a Wronskian, the 

eigenvalues with ^(z) > are precisely the roots of 822(2). By similar arguments, the eigenvalues in the 
open lower half-plane are precisely the roots of 8ii(z). There are no real eigenvalues, because according to 
Proposition A.6 all solutions oscillate for large |a;| when z is real. Suppose that z is an eigenvalue in the 



upper half-plane. Then it follows that there is a nonzero proportionality constant rj such that 
(180) ii(x;z) ^r]i+(x;z). 

Let /, //, ///, and IV be the four quadrants of the plane. The following corollary can be obtained 



from Proposition A. 8 since 5*11(2) extends from the real axis to the lower half-plane as 822(2* 



Corollary A. 9. // z is an eigenvalue on the imaginary axis, then —z is also an eigenvalue. Similarly, if 
z € I is an eigenvalue, then —z* G //, — z G ///, and z* G IV are also eigenvalues. 

As a result, the eigenvalues come either in pairs (on the imaginary axis) or in quartets off the axes. Using 
the symmetries 



(181) 



ii(x; z) 



1 
-1 



j^(a;;-z) = 



1 

-1 



J2 (x;z ) 



(which follow from ( 173 1 and ( 178 ) upon extension to complex z) the relation ( |180 1 holding for an eigenvalue 
z with 3(z) > implies that 



(182) 



ji (a;;-z*) = ?7*j2 (a;;-z*), (x; -z) 



-vit(x; 



-Z), i2(^:2*) ^ -V*il(X:Z*)■ 



l^ z is an eigenvalue on the positive imaginary axis, then these symmetries show G M. Note that if 82i(z) 
has a meromorphic extension into the upper half-plane and is finite and nonzero at an eigenvalue z in the 
upper half-plane, then 77 — 82i(z). 

Definition A. 10. Suppose that 822(z) has only simple zeros in the open upper half-plane. The scattering 
data for the Cauchy problem consist of (i) the reflection coefficient 

82i(z) 



(183) 



P(2) 



z G 



822(z) ' 

(ii) the eigenvalues, or the zeros {z„} of 822(z) in the open upper half-plane, and (Hi) the modified propor- 
tionality constants {c°} where 

Vn 



(184) 



822(Zn) 



It turns out that this information is sufficient to reconstruct the potentials / and g, assuming that ^22(2) 
has no real zeros or complex multiple zeros. 

Up to this point, we have considered (1231 (or equivalently (127l) for u — f(x) and eut — g(x) as 



fixed functions of x. However, if u(x,t) evolves in time t according to the sine-Gordon equation, then / 
and g will depend parametrically on t, and so will the Jost solutions of ( |127[ ). We must therefore expect 
the scattering matrix S(z) = S(z;t) and the proportionality constants — ijn(t)} to vary with t. Since 
the sine-Gordon equation is the compatibility condition for (1231 and ( |124 1, we may use (1241 to calculate 
the time dependence. In order that all quantities of interest remain well-defined as t varies, we introduce a 
technical condition on solutions u(x,t) of the sine-Gordon equation ([T]). 

Definition A. 11. Let p >1. A solution u(x,t) of the sine-Gordon equation ([T|) is said to have U'-Soholev 
regularity if I — cos(u) sin(u), u^, u^x, Ut, md Utx all exist in the sense of distributions and lie in the space 
L^(IR) as functions of x for all t > 0. 
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Appendix [B] contains a proof of the fact that as a dynamical system, the sine-Gordon equation 
preserves L^-Sobolev regularity, so in fact, it is only a condition on initial data. The case of interest in 
inverse-scattering theory is p = 1 . 

Proposition A. 12. Suppose that u — u{x,t) is a solution of the sine-Gordon equation ([T|) having L^- 
Sobolev regularity. Then the corresponding time evolution of the scattering data computed by solving (1271 
with potentials f{x) = u(x,t) and g{x) = eut{x,t) is given by 



(185) 
(186) 
(187) 



iS'ii(z) and S22{z) (and thus the eigenvalues) are independent oft, 
Si2iz; t) = 5i2(z; 0)e-2'^(-)*/^ S2i{z; t) = S2i{z- 0)e2*^(-)*/% 
Vn{t) = 77„(0)e2'^(-")*/^ 



where rjn is the proportionality constant associated to the eigenvalue z„ in the open upper half-plane. 
Proof. Since u{x,t) is a solution of the sine-Gordon equation, we can find functions Ci{t;z), cj(i;z) (inde- 



pendent of x) such that :— c^ji and cjj2 simultaneous solutions of the Lax pair (1231 and 

( 124[). Inserting into (124) and using the relation (1301 between jj^ and m^, we find 



(188) 



4je— log{ci) m.^ + Aiem.'^^ = Bni]^. 



The limits of and B as a; ^ -l-oo exist by Proposition A. 6 and equation (1551. We now show m^^ also 
has a limit as x ^ -l-oo. Taking a time derivative of ( 149 1 shows 

(189) m+{x;z,t) = - K^^{y■,z,t)m+{y■,z,t)dy - {y; z,t)m+{y; z,t) dy . 

J X J X 

For z e {3(z) < 0} n {\z\ > 6} for any fixed 5, we have Kj'^, K.^^ E and mf E as functions of y by 
the assumptions that 1 — cos{u),sm{u),Ux,Ut,UxxTUtx G (note uu = u^x — sin(M)/£^). Therefore by an 
iteration argument, m^^ E L°° as a function of x uniformly for z E {'^{z) < 0} n {\z\ > d}. By an analogous 
argument using the zero gauge, one sees that mfj E L°° as a function of x for 3(z) < 0. It then follows from 
a dominated convergence argument applied to (1891 with this new information in hand that 



(190) 



lim m^((a;; z, t) — 0. 

x — ^+oo 

a(2)<o 



-cG gives 



"1' 






= 4D 









Using this result to take the limit of ( 188 1 as a; 

(191) 4ie^log(c+) 

and so up to a multiplicative constant (independent of x and i), "w^ — e^*^*/^j]^ for 5(z) < 0. Similar 
arguments for the other Jost functions show that in the respective closed half-planes of existence 

(192) w± = e-'-°*/^jf and w± = e'-°*/^j^ , 



again up to multiplicative constants. Substituting the expressions ( 192 1 into the time-evolution equation 
( fl24t gives in particular 



(193) 



ejf = eBJ^ + iDJ^a^, . 



z E 



Differentiating equation (167) gives 
(194) 



= Jt S + J St. 



Substituting (193) into (|194| and using (|167|) gives 
(195) 



£-S(z) - zD{z)[S{z),a3] = 2D{z){S{z)a3 - a^Siz)) 



which yields (185) and (186) 
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For the time evolution of a proportionality constant rj associated to an eigenvalue z in the upper 



half-plane via (180), we differentiate the relation (|180| with respect to t and obtain 
(196) 



ht(^;z,t) = J2 {x;z,t}— +?7J2t(a;;z,t) 



Obtaining the time evolution of {x; z,t) by substituting (1921 into ( 124[ ) gives 

(197) eBji {x;z,t) + iD{z)i-i^ {x;z,t) === ej^(x; z, i) ^ + er?Bj^(x; z, i) - iriDiz)^^ {x; z,t) 
Using ( 180 1 to eliminate j;^ {x; z, t) gives, since j^(x; z, t) ^ [0, 0]^, 

(198) = 2iD{z)ij, 

at 



which gives equation (1871. 



□ 



A. 3. The matrix G and its properties. We now introduce a piecewise-meromorphic function G{z) = 
G(z; x,t,e) whose singularities and jump discontinuities encode the scattering data. Assuming that u is a 
solution of the sine-Gordon equation with L^-Sobolev regularity, define a corresponding matrix G(z) by 
(199) 



G(z) := 



-5*22(2) 



m.^ (x; z, t), m^(x; z, i) 



S22{z) 



Ji (x;z,i)e^^-/M+(x;z,t)e-^^-/^ 



m+(x;z,i), 



S'ii(z) 



m2 ix;z,t) 



itix; z, ^)e'^-/^ -i— j2 {x; z, t)e-^^^/' 



3(z) > 
3(z) < 0. 



From the symmetries of the scattering matrix (Proposition A. 8 1 and of the Jost functions (181 1, G(z) has 
the symmetries 

(200) G(z) = tT2G(-z)a2, G(z)-i = -G(z*)t. 

The matrix G(z) will have poles at the eigenvalues {z„} due to the presence of Sjj{z) in the denominator. 

Proposition A. 13. Suppose that u is a solution of the sine-Gordon equation with L^-Sobolev regularity, and 
let G(z) be the corresponding matrix defined by (199 1. Let Zn £ I be a simple eigenvalue (that is, a simple 
root of S22{z)) and let rjn be the corresponding proportionality constant defined by equation {180). Then 

0" 



(201) 

where 
(202) 



ResG(z) = lim G(z) 
ResG(z)= lim G(z) 



a 

c„ 




ResG(z)= lim G(z) 



ResG(z) = lim G(z) 








„ _ J) p2iE{z^)x/e+2iD{z^)t/E „0 



^n(O) 
S22i.Zn) 



These formulae also hold for eigenvalues Zn on the positive imaginary axis, in which case we have c„ S iM. 

first quad 

ResG(z) 



Proof. Consider z„ in the first quadrant and let O2 := [0,0]^. Using (180 1 

1 



Ji (x;z„,^)e^^(^'•)-/^02 



(203) 



^22^^-0) 
r?„(t)e2*-^(^")^/= 



^22^'^ri) 



Jt(x;z„,^)e-'^(^"W^02 



lim G(z) 





Cn 



The other three proportionality constants are handled similarly with the help of (182). 



□ 



The matrix G(z) also generally has a jump discontinuity across the real axis. This discontinuity 
is expressed in terms of the reflection coefficient defined by (183). Note that, by definition, p(z) = 
for refiectionless potentials. Note also that if p(z; 0) = 0, then p{z] t) = for all i > according to 
Proposition |A.12[ 
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Proposition A. 14 (Cheng et al. [1]). Suppose that u is a solution of the sine-Gordon equation with L^- 
Sobolev regularity, and let G{z) be the corresponding matrix defined for z e C \ M by (1991. Suppose also 
that Sii{z) and S22{z) do not vanish for any real z. Then Giz) takes continuous boundary values G+{z) 
and G_(z) for z G M from the upper and lower half-planes respectively. Moreover, these boundary values are 
related by G+(z) = G_(z)V(z), where the jump matrix is defined for z € M as the continuous function 

1 + \p{z; 0) |2 -p{z- o)*e-^<E^+Dt)/e- 



V(z) ^Y{z;x,t,s) := 



(204) 



The jump matrix satisfies 
(205) 



1 -p(z;0)*e-2^(^^+^*)/^' 
1 



lim V(z) = lim V(z) 



1 0' 

-p(z;0)e2'(^"+^*)/^ 1 



Proof. The continuity of the boundary values follows from the definition (199) and established properties of 
the columns of M , given that the denominators Sjj{z) do not vanish for z G M. 
From equation (167l and det(S(z)) = 1, 

(206) 











'5'22 


— 'S'12 


Jl 7J2 




Jl 1J2 




'S'21 


Sii 



Therefore, using (1301, we obtain 
1 ^- 



(207) 

and 
(208) 



1 

-^"12 = -^J2 e 

»Jii "Jii 



-lEx/e ^ I- c 1+ I c '■ + \^-iExle '^12 „-2iEx/e 
— (,-'J12Ji + J11J2 je — p 



Sii 



m. 



^11 



1 ^- 



'522 '-'22 



■jTe^^-/^ = ^(^22jt - 52ut)e^^^/^ = m+ - ^e^^^^/^mt- 

'J22 '^22 



Combining these two equations in matrix form using the identities 



(209) 
gives 

(210) 



G+(z) [S'22(2:) ^m.^{x;z,t),Ya^{x;z,t)], G_(z) = [m]^(a;; z, t), 5ii(z) ^Ya.^{x;z,t)] 



G+(z) 



1 



p{z;t) 



,2iEx/e 



G_(z) 



1 -p(z;t)*e-2«-^"^/=' 
1 



which, after taking into account P ropo sition A. 12 shows (2041 
Finally, using Proposition 



A.6 



the Wronskian relations ( 168 1 shows that 
(211) 



the relations (|130 1 between the columns of M and those of J , and 
[, lim S(z) = I, 



limS(z) = (-l) 

2^0 



which implies the limits (205) of the jump matrix. 



□ 



Therefore in any reflectionless case the condition p = implies that V(z) = 0, and by continuity 
of the boundary values holding in the absence of real zeros of Sjj{z), a Cauchy integral argument shows 
that G(z) is meromorphic in the entire z-plane. The next proposition follows immediately from (159) in 
Proposition A.6 and (211 1. 

Proposition A. 15. Suppose that u is a solution of the sine-Gordon equation with L^-Sobolev regularity. 
Then the matrix G(z) defined by (199) satisfies 

(212) lim G(z) = I . 

Z — >OQ 

(The limit is independent of direction in the complex plane.) 

The next proposition describes how to recover cos(u) and sin(u) (and thus u) from G(z). 
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Proposition A. 16. Suppose that u is a solution of the sine-Gordon equation with L^-Soholev regularity. 

.(0) 



Let G'-'^' :— lim G{z). Then cos(m) and sin(u) are expressed in terms of the elements G^^ of G^'^^ by 



(213) 

and 

(214) 



cos(u) 



^(0)^(0) 



^(0)^(0) _ , . ry^yu) ^y\J 

<-ri2 ^21 — + ^<-^12 ^21 



(0)^(0) 



sin(u) =2G^"^G^2^ 



9^(0)^^(0) 



Proof. From Proposition A. 14 we have V(0) = I, and it follows that G^'^^ is well defined. Also, from Lemma 
A. 5 Gx{z) remains bounded in the limit z 0. Since the columns of G{z) satisfy the differential equation 
( 127[ ) where they are defined, 

(215) AieG^ = iEcTsG - iEGa^ + QG. 

Let 



(216) 



Q 



(-1) 



lim zQ 



- cos(m) 

- sin(u) 



— sin(ii) 
-(1 — cos(u)) 



Then multiplying (215 1 by z and taking the limit as z — > gives 

(217) 0= -a3G(°) + G(0V3 + Q^"'^G(°). 

Therefore 

^(-1) = ^3 _ G(")a3G(o)-\ 



(218) 



Q 



which gives (2131 and (214). The consistency and reality follow from the holomorphic and antiholomorphic 
symmetries (2001 at z = 0. The Pythagorean identity follows from the identity det G(2) = 1, which in turn 
follows from (|130| and ([T68|. □ 



A. 4. Riemann-Hilbert problem. From this point onward, we assume that u = u{x, t) is a solution of the 
sine- Gordon equation with the following properties: 

(a) The solution u has L^-Sobolev regularity. 

(b) All eigenvalues corresponding to u are simple. 

(c) The entries 5*11(2) and S22{z) of the scattering matrix S(z) have no real zeros. 
These are natural conditions|fl 



While G(z) has been defined (see (199l) in terms of solutions of differential equations with non- 



constant coefficients involving u(x,t), we have shown directly that G(z) has certain properties that depend 
only on the scattering data for u, which may be calculated at any convenient instant of time, say, t = Q. 
The Riemann-Hilbert problem asks whether it is in fact possible to determine G(z) (and hence u, via Propo- 

To pose the Riemann-Hilbert problem, we seek a matrix 



sit ion 
H(z) 



A. 16) purely from the scattering data. 



(219) 



H(z; X, t, e) with the following properties: 
(I) H(z) — > I as z ^ 00 uniformly in all directions. 

(II) H(z) is meromorphic in the upper and lower half-planes with simple poles only. The residues of 

H(z) are required to satisfy (201 1 and (202). 
(Ill) On the real axis z S M, H(z) obtains continuous boundary values and satisfies the multiplicative 
jump relation 

H+(z) = H_(z)V(z) 



with V(z) given by ( |204 l. 

Note that the problem description only involves the scattering data, and also that x and t enter in a simple 
explicit manner. 



^According to the arguments to be given in Appendix [b] the solution to the Cauchy problem corresponding to the special 
initial data |2jl satisfies condition (a) for all S M. Note, however, that for fixed fi ^ 0, condition (b) does not hold for 
this data if e satisfies {\/^^ + 1 — l)/2s G {1,2,3, . . . }. Furthermore, for any fixed fi, condition (c) fails for this data if e 
satisfies (-^/x^ -I- 1 — l)/2e S {1, 2, 3, . . . } as e J, 0. In other words, conditions (b) and (c) fail infinitely often in the limit e J, 0. 
Nevertheless, conditions (b) and (c) both hold generically with respect to e. Also, if conditions (b) and (c) hold at a fixed time 
(say t = 0), they will hold for all time. 
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Theorem A. 17. Assume conditions (a), (b), and (c). Then the Riemann-Hilbert problem has a unique 
solution for all x Cz and t > 0. 



Proof. As a consequence of the L^-Sobolev regularity of the solution to the Cauchy problem, the matrix G{z) 
defined by the formula (199) can in principle be constructed for any a; e M and t > 0. From Propositions 



A. 13 A. 14 and A. 15 G{z) is, therefore, a particular solution of the Riemann-Hilbert problem. 

It remains to determine whether the solution is unique. First, suppose that H(z) is any solution to 
the Riemann-Hilbert problem. We will show that det(H(z)) = 1. Indeed, suppose zq is a (simple) pole of 
H(z) with associated residue parameter cq. From the residue condition (II), H(z) has a Laurent series of 
the form 



(220) 



H(z) 



z - Zo 



Y.hi-\z^zoy 



n=0 



where h^" are constant vectors. Expanding the determinant by columns shows det(H(z)) — 0{1) near 
z — Zq, that is, det(H(z)) has no poles, and so for all z ^ M det(H(2:)) is a scalar analytic function 
of z. Moreover, from the jump condition (II), det(H_|_(z)) = det(H_(z)), and also det(H(z)) obtains 
continuous boundary values along the real axis. Since det(H(z)) is analytic in the open upper and lower 
half-planes, Cauchy's theorem shows det(H(2)) is analytic on the real axis as well. Therefore det(H(z)) is 
an entire function, and from the normalization condition (HI) we have det(H(z)) —^ I as z —^ oo. Therefore 
det(H(z)) = 1 by Liouville's theorem. 

Now assume H(z) is another solution to the Riemann-Hilbert problem and consider the matrix 
E(z) :— H(z)H(z)^^. Using det(H(z)) = 1, direct multiplication shows that all singularities of E(z) are 
removable. For z e M, 

(221) E+(z) = H+(z)H+(z)-i = H_(z)V(z)V(z)-iH_(z)-i = H_(z)H_(z)-i = E_(z), 



so E(z) has no jump discontinuity across the real axis. Furthermore, E(2) also achieves continuous boundary 
values on the real axis, and therefore, by a Cauchy argument, is an entire function. In addition, E(z) I 
as z ^ oo as this is true of both H(z) and H(z). Therefore, up to removable singularities, E(z) is entire 
and bounded in the complex plane, and again by Liouville's theorem E(z) = I, or put another way, H(z) = 
H(z). □ 



The significance of this result is that the solution of the Cauchy problem for the sine-Gordon equation 
corresponding to a broad class of initial data can be completely characterized through the solution of the 
Riemann-Hilbert problem given the scattering data calculated at i = (which therefore encode the Cauchy 
data / and g). In particular, this point of view is well-suited to asymptotic analysis in various limits of 
interest, for example, e J, 0. 

Note that, in formulation presented above, the existence of a classical solution H(z) = G(z) to 
the Riemann-Hilbert problem follows from the well-posedness theory for the sine-Gordon Cauchy problem 
explained in Appendix |B] In the absence of such an independently developed theory of well-posedness, it is 
still possible to prove the existence of a solution to the Riemann-Hilbert problem by direct means in various 
situations. See, for example, [27] . 

To close this appendix, we comment on the semiclassical asymptotic analysis of the Riemann-Hilbert 
problem corresponding to the special initial data (|2]), which is work in progress. The so-called steepest descent 
method for matrix-valued Riemann-Hilbert problems developed by Deift, Zhou, and their coworkers is a pow- 
erful method of asymptotic analysis in which a basic technique is to introduce a sequence of carefully-chosen 
piecewise meromorphic changes of the dependent variable H(z) (frequently these are called deformations of 
a Riemann-Hilbert problem). The ultimate aim of these deformations is to move the jump discontinuities 
from one contour to another in order to exchange oscillations for exponential decay. A simple example of 
such a deformation arises from the factorization (204) of the jump matrix V(z) (see, for example, [25] ). 
The idea is to replace H(z) with another unknown H^^\z) defined as follows. Let D± be bounded subsets 
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of C± both adjacent to the same interval / of the real axis. Then we set 



H(i)(z) := H(z) 

(222) 

H(i)(z) H(z) 



1 0' 

pl^^.Q-^^2t{Ex+Dt)/e I 

1 -p(^*;0)*e-2«(i?-+ct)/e- 
1 



zeD+, 



and, for all other z in the upper and lower half-planes, we simply set H(^)(z) H(z). Note, in particular, 
that as a result of the factorization (2041, the matrix Il^^\z) extends continuously to the real interval /, that 



is, there is no longer any jump discontinuity across /. Typically, the residues of the poles of H(z) = G{z) 
are exactly cancelled in the regions D± by this deformation. That is, a deformation based on ( 204[ ) has the 



added benefit of removing the pole^from the Riemann-Hilbert problem. However, if S2i{z) has (phantom) 
poles (as in the case corresponding to the initial data (l2|; see the discussion just before the statement of 



Proposition 2.9 1 and if any of these poles lie in the domains D±, then the above deformation will also 
introduce new poles into the matrix H(^^(z) where there were none in H(z). Avoiding the phantom poles 
places additional conditions on the regions D± in which the change of variables can be made, and these 
conditions further complicate the steepest-descent analysis. 

Appendix B. The L^'-Sobolev theory of the Cauchy problem for the sine-Gordon equation 
Consider the Cauchy problem corresponding to topological charge n: 

{PDE: e^utt- e'^Uxx+sm{u) ^0, xeR, t>0 

Initial Conditions: u{x,0) = f{x) , sut{x,0) = g{x) 

Boundary Conditions: u{—oo, t) = , u{+qo, t) = 27rn , n £ Z . 

Here /(•) and g{-) are given functions; / satisfies the given boundary conditions and g vanishes as |a;| oo 
in a sense to be prescribed later. Regardless of the value of the topological charge, we may easily transform 
this problem into one for which the corresponding boundary conditions are homogeneous by making the 
substitutions 

(224) u(x, t) = U{x, t) + b{x) , ,f{x) = F{x) + b{x) , 

where &(•) is function that satisfies the (typically nonhomogeneous) boundary conditions. For convenience, 
we will suppose that b'{-) is in the Schwartz space, and moreover, that h{—x) — 27m = —b{~x). A concrete 
example of a function 6 of this type is given by 

(225) h{x) := 7rn(l -|- tanh(a;)) , 

but our arguments will not rely on this particular definition. Thus, U and F satisfy homogeneous boundary 
conditions, and so we arrive at the equivalent Cauchy problem: 

fPDE: Utt-Uxx = Q{x,t), xeR, t>0 

(226) < Initial Conditions: U{x, 0) = F{x) , Ut{x, 0) = G{x) 
[Boundary Conditions: C/(— oo,i) = 0, U{+oo,t) = 0, 

where 

(227) Q{x, t) := b"{x) - \ sm{U{x, t) + b{x)) 
and, for symmetry of notation, we have defined 

(228) G{x) := -^g{x) . 

^Of course, in the (nongeneric) reflectionless cases for which e = EjvCm) hence p = making 'V{z) = I, both factors of 
V{z) as written in | |204| are trivial, so the poles may not be removed in this way. However, the poles may indeed be removed 
by finding an appropriate analytic interpolant of the proportionality constants at the corresponding eigenvalues {zn} and 

using this interpolant along with the formula ( |104| to construct a meromorphic function whose residues at the eigenvalues 
are the modified proportionality constants {cn}- In I12| . this approach was used to remove the poles from the refiectionless 
inverse-scattering problem associated with Klaus-Shaw initial data for the focusing NLS equation. 
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We may express U in terms of the initial data F, G, and the source term Q with the help of Green's 



function: 



(229) u{x,t)^^Fix + t) + ^F{x-t) + ^ 



x+t 

x~t ^ Jo 



t fX + it-to) 



X-(t-to) 



Q{xo,ta)dxodto. 



Since Q depends on U, this formula (Duhamel's formula) does not represent the solution in closed form, but 
as derivatives of U are not involved it amounts to a useful reformulation of the Cauchy problem ( 226 1 in 
weak form. We introduce a nonlinear operator T by 



(230) 
where 



T[U]{x,t) -.^ A{x,t) + S[U]{x,t) , 



(231) A{x,t):^^Fix + t) + ^Fix-t) + ^ 



x+t 



G{xo) dxo 



t fX + it-to) 



b" (xq) dxo dto 



Jx-(t-to) 



and 
(232) 



S[U]{x,t) ■.= -—^ 



t fx+{t-ta) 



Jx-(t-tn) 



sin([/(xo, to) + 6(2:0)) dxo dto ■ 



Thus, the weak formulation (2291 of the Cauchy problem (226) takes the form of a fixed-point equation in 
function space: U — T[U]. 

We want to study this problem in the spaces with respect to x. For each T > and p > 1, define 
a norm by 



(233) 



\U\\xS, ■= sup \\U{-,t)\\LP^ sup 

0<t<T 0<t<T 



+00 



\U{x,t)\Pdx 



and a corresponding Banach space of functions U : 

(234) Xf, := |f/(x,t), xeR,t£ [0,T] such that \\U\\x^ < ooj 



Lemma B.l. Suppose that F e LP(M.) and G G LP{M). Then the function A defined by (231) lies in the 
space for each T > 0, and 



(235) 



\A\\x-, < \\F\\l.+T\\G\\l. + '^W\\l. 



Proof. By the Minkowski inequality and translation invariance of L^, 

\\A{;t)U.<\\FU. + ^ 

(236) 



4- 00 



x+t 



G{xo) dxo 



dx 



+ 



1 



+00 



t i-x+(t-to) 



b"{xo) dxo dto 



Jx~{t-to) 



dx 



The integral triangle inequality gives 



(237) 



ll^(-,t)l|L. <I1^I1l.+ 
1 

+ 2 



-00 / px+t \ P 

|G(a;o)| dxo I dx 



00 \J x—t 



+00 / j-t px+(t-to) \P 

/ / \b"{xo)\dxo dto dx 

00 Wo Jx-{t~to) I 
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Now, by Holder's inequality we have 

(■X+t 

[[ 

(238) 



and 
(239) 





■ px+t 


1 

Q 


H 


/ 1« dxo 






Jx-t 





rx-i-i 

/ \GixoW 

Jx-t 

= (2t)f / ^ \Gixo)\Pdxo 

Jx-t 



dxo 



p 



/ \b"{xo)\dxodto] < 

Jx-(t-ta) I 



t fX+it-to) 



1^ dxo dto 



Jo Jx 

where q satisfies l/p+l/q = 1. Therefore, 



Jx-{t-to) 
2p (-1. (-x+it-to) 



t fX + it-tf,) 



11 

Jo Jx 



\b"{xoWdxodto 



Jx-(t-to) 



1\ P 
P 



X—{t — to) 



\b"{xo)fdxodto 



(240) 



\\A{.,t)\\L^<\\FU. + l{2t)l 



/ \G{xo)\Pdxodx 

Jx-t 

b"{xo)\P dxQ dto dx 



1 2 
+ 2*^ 



oo Jx-t 

+ X fl fX + it-to) 



-OO Jo J X- 



[t-to) 



Using Pubini's theorem to exchange the order of integration, we have 



OC fX-\-t 



|G(xo)|^ dxQ dx 



DC fXi:)+t 



(241) 



and 



(242) 



OO J XQ—t 
+ 00 



/i-oo rXQ-\-z 
\G{xo)\P / dx 
-OO J Xn—t 



\G{xo)\P dx dxo 

Xo+t 

xo—t 



dxo 



nG\\%, 



+00 i-x+(t-to) 



/i-oo j-l. r- 
-OO Jo Jx 



-{t-to) 



\b"{xo)\Pdxodtodx = 



t Z' + OC />X-\-{t — to) 



11 

pt n-\-oo p 
Jo J -OO Jx 



-OO J x—{t — to) 
f^t p-\-oo pxo-\-{t—to) 



\b"{xo)\Pdxodxdto 



xo — {t—to) 



\b"{xo)\Pdxdxodto 



LP ■ 



\\A{;t)\\LP < \\F\\l. +t\\G\\LP + ^nb"\\L. , 



Therefore, using l/p+ 1/q = 1, we have 
(243) 

and taking a supremum over t G [0, T] completes the proof. □ 

In terms of the function b{-), let another function c(-) be defined as follows: 

c{x) := \b{x)\x-{x) + \b{-x)\x+{x) , 

where x+i^) and X-i^) are the characteristic (indicator) functions of the sets a; > and a; < 0, respectively. 
This function is bounded and rapidly decaying as |a;| — > oo. 

Lemma B.2. Let T > 0, and suppose that U G X^. Then S[U] G X^, and 



(244) 
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Proof. By the triangle inequality for integrals we have 



\\S[U]{.,t)\\L^ = 



(245) 



< 



2e2 

1 
2? 



+00 



+00 



sin(J7(a;o,to) + b{^o)) dxo dto 



Jx-{t-t„) 

/, fX + {t-to) 



dx 



I sin(L''(a;o,io) + ^(^o))| dxo dto 1 dx 

lO Jx-{t-to) ) 

For any real x and t we have, using periodicity of the sine function, the inequality | sin(a;)| < |a;|, and the 
property b{x) — 2nn = —b{—x) that 

|sin(;7(a;, t) + b{x))\ = \sm{U{x,t) + b{x))x-{x) + sm{U{x,t) - b{-x))x+{x)\ 

< \sm{U{x,t) + b{x))\ x-{x) + \sm{U{x,t) - b{-x))\x+{x) 

(246) < I U{x, t) + b{x) I X- {x) + I U{x, t) - b{~x) \ x+ {x) 

<\U{x,t)\ + \b{x)\x-{x) + \b{-x)\x+{x) 
= \U{x,t)\+c{x). 

Therefore, by the Minkowski inequality. 



\S[U]{;t)\\L. < 



(247) 



1 

2^2 



/ \U{xo,to)\dxodto dx 

Jx-(t-to) I 



1 

2? 



+00 / |.t r.X + {t-to) \^ 

/ / c{xq) dxo dto I dx 

-oo \Jo Jx-{t-to) I 



Applying Holder's inequality, we then find 



\\S[U]{ 



< 



(248) 



2 

t1 

2i2 



+ 



+ 0C fl. fX + (t-to) 



11 

J —CO JO Jx 



2 
tl 
2i2 



X— (t — to) 
+00 i^t rX-\-{t—to) 



I U {xq , to) r dto dx 



-oo Jo J X- 



(t-to) 



c{xo)^ dxo dto dx 



where 1/p -\- 1/q = 1. The order of integration may be exchanged by Fubini's theorem: 



+00 pt pX-\-{t—to) 



/-\-oo pt p 
-00 Jo Jx 



-{t-to) 



\U{xQ^tQ)\^ dxQ dto dx ■ 



t p-\-oo pX-\-{t—to) 



pz p-\-oo p 
Jo J — 00 Jx 



(249) 



X-{t-to) 

t f+ca i-xo + it—to) 



\U{xo, to)\^ dxo dx dto 



J — 00 J xo — {t—to) 
t 



\U{xo, to)Y' dx dxo dto 



Similarly, 
(250) 



= / ||?7(-,io)||ip-2(i-to)rfio 
Jo 

Kt" sup \\U{;to)\\l.. 
0<to<t 

/ + OO pt rx + {t — to) 
I I c{xoY dxodtodx = t^^cYLP ■ 

-00 Jo Jx-{t-to) 



The proof is therefore complete upon taking a supremum over < t < T. 
Lemma B.3. Let T > 0. Then whenever U and V 



(251) 



mU]-T[V]U.<^^\\U-VU.^. 



□ 
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Proof. Clearly, we have T[U] — T[V] = S[U] — S[V]. (We could have assumed further that F and G lie in 
L^, so that T[U] and T[V] are individually well-defined as elements of X^^, but as A{x, t) cancels out of the 
difference, it is not necessary to make such an assumption here.) Now by the triangle inequality for integrals. 



\\Sm-,t)~S[V]{;t)\\L. 



1 

2? 



+00 



:+(i-to) 



{sin(C/(a;o,to) + b{xo)) 



Jx~{t-to) 



sin(y(a;o, to) + ^(xo))} dx^ dto 



(252) 



dx 



< 



2e2 



+00 



Jx-{t-to) 



|sin(C/(xo,to) + b{xo)) 



— sm(y{xQ,tQ) +b{xQ))\ dxQdtQj dx 
Now, since by the mean value theorem, | sin(a;) — sin(j/)| = | cos(^)| • |a; — y| < |a; — y|, we have simply 



(253) \\S[U]{;t)-S[V]{;t)\\L^< 



2e2 



+00 



Jx-(t-to) 



\U{xa,tQ) -V{xQ,ta)\dxQdtQ\ dx 



The proof is then finished upon using Holder's inequality and Fubini's theorem in exactly the same way as 
in the proofs of Lemma |B.1| and Lemma |B.2| □ 

Theorem B.4 (Local Existence and Uniqueness in L^). Let p > \ and suppose that F,G ^ L^. Define 
T := min{e, 1}. Then there exists a unique weak solution of the sine-Gordon Cauchy problem ( 226[ ) (that is 
a unique solution of (229 I J in the space XJ^. 



Proof. Combining the results of Lemma B.l and Lemma B.2 we have that for any U € X, 



(254) 

Since T < 1, it is also true that 



nmwx- < \\f\\lp + \\g\\l.t+ ^iife'iiLpT' 



\nu]\\ 



< \\F\ 



1, 



LP 



1 

2? 



1 



\c\\lpT' 



rp2 

2? 



2e2 



\F\ 



LP 



\G\\ 



1, 



LP 



(255) 
Let 

(256) R = Rp[F,G]:^2 
and let Bfi denote the bounded subset of Xi^ given by 

(257) Br:= e X^ such that \\U\\xp < i?} 



y2 

2? 



m 



m 



LP 



1 

2? 



C\ LP 



Then, (2551 takes the form 
(258) 



< 



R T2 



Mx'^ ^ R 

2 - 2 



\U\\x^^ 



2 e2 2 - 2 2 ' 

where we have used T < e. Therefore, U S Bj^ implies that T[U] G Bji as well; that is, T maps the bounded 
set Bn to itself. Furthermore, combining the inequality T < e with Lemma |B.3| we obtain 



(259) 



mU]^T[V]\\j,P < -\\U-V\\j,P, 



U,VeBR. 



(Actually, this holds for all U,V E Xlj^.) Therefore, T defines a contraction mapping on Br, and so there 
exists a unique fixed point U G Br of the mapping T, that is, a unique solution in Br of the equation 
U — T[U] equivalent to the weak form (2291 of the Cauchy problem. It is easy to see that the number R 
could also have been replaced by any larger number, and therefore the solution guaranteed by the contraction 
mapping principle is actually unique in the whole space X^. □ 
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Theorem B.5 (Global Existence and Uniqueness in L^). Let p > I and suppose that F, G, and F' all lie 
in LP. Then for each finite T > there exists a unique weak solution of the sine-Gordon Gauchy problem 
(226 1 in the space X^^. 

Proof. We wish to iterate the argument in the proof of Theorem |B.4| by restarting the Cauchy problem at 
time T with new initial data for which F{x) is replaced by Fi(x) := U{x,T) and G{x) is replaced by the 
distributional derivative 6*1(2;) := Ut{x,T). Since 



(260) 



If^Wl. - \\u{;T)\\lp < ||;7|Up <oo, 



the new initial data satisfies Fi E L'p. However, the norm of does not directly provide us with any control 
of ^-derivatives needed to place Gi in the space L^ together with Fi. 
To analyze Gi, differentiate (2291 with respect to t: 



(261) 



Ut{x, t) = ^F'{x + t)- ^F'{x -t) + ^G{x + t) + ^G{x - t) 
+ \j [b"{x + {t- to)) + b"{x -{t- to))] dto 



[sin{U{x + (i - to), to) + b{x +{t~ to))) 

+ sm{U{x- {t-to),to) + bix~ (t-to)))] dto . 



Jo 



It follows that 



\\Ut{;t)\\L. < \\F'\\lp + \\G\\l. 



(262) 



1 
2 
1 
2 

1 
2? 

1 
2? 



b"{- + {t-ta))dto 



b"{--{t-to))dto 



LP 



LP 



sin(C/(- + {t~ to), to) + hi- + it- to))) dto 



sin(C/(- - (t - to), to) + b{- - (t - to))) dto 



LP 



LP 



Now, by Holder and Fubini arguments, 



b"{-±{t^to))dto 



LP 



(263) 



< 



1 

< ti 



+ C30 



+ C30 

-00 \J0 

-00 pt 



00 Jo 
t /'4-00 



1 

tl I ^ 

J-00 

tP^h\b"\\LP 
t\\b"\\Lp. 
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b" {x ± {t - to)) dto dx 

\b"{x±{t-to))\dto] dx 
b"{x±{t-to))\Pdtodx 
b"{x±{t-to))\Pdxdto 



Using as well ( 246 ) and the Minkowski inequality, 



(264) 



sin(C/(- ± (i - tn),to) + b{- ± {t - to))) dto 







LP 



< 



< 



< 



1 

< t1 



1 

= t1 



+ 00 



+ 00 



+ CX2 



sin{U{x ± (i - to), io) + ± {t - to))) dto 



dx 



\sm{U{x±{t-to),to) + b{x±{t-to)))\dto dx 



\U{x±{t-to),to)\dto+ c{x±{t-to))dto] dx 



\Uix±{t-to),to)\dto] dx 



+ 00 



oo Jo 



\U{x±{t-to),to)\^ dto dx 



+00 



\U{x ± {t ~ to),to)\P dx dto 



1 

p 1 
+ t'} 

1 

p 1 
+ ti 



+ 00 

-oo Jo 

t /' + 00 



c{x ± {t — to)) dto j dx 

c{x ± {t - to))^ dto dx 
c{x ± {t - to))^ dx dto 



<t sup \\U{-,to)\\LP+t\\c\\LP 
0<to<t 



Therefore, 
(265) 



\\Ut{;t)\\LP < WF'Wlp + \\G\\lp + t\\b"\\LP + ^\\c\\lp + ^ sup \\U{-,to)\\LP 



Thus, if the initial data (already assumed to satisfy F ^ and G G to guarantee the existence oi U G 
for t G [0,T] according to Theorem B.4| also satisfy F' e L^, then [/^ is uniformly in for all i e [0,T], 
and in particular the new initial condition satisfies Gi £ . This is sufficient to iterate the argument in the 
proof of Theorem |B . 4| an arbitrary number of times, with a fixed time step T, and the proof is complete. □ 



Theorem B.6 (Global L^-Sobolev Regularity). Suppose the same conditions as in Theorem B.5 namely 
that F, F' , and G are all in L^iK). Then the unique global weak solution of the Cauchy problem (226 1 
satisfies U € LJ^JLP(M)), € L^^{LP(R)), and Ut G L^,{,LP{ 



Moreover, if the initial data have one 
more derivative in (that is, F" and G' are in as well as F, F' , and G), then this further regularity is 
preserved as well: one also has Uxx S -^i^c(^^('^)) '^'^'^ ^ -^iSc(^^('^))- 



Proof. The fact that U G Lj^^(LP(M)) follows from the statement of Theorem B.5 and the fact that Ut G 
LJ^^(L^(R)) follows from the estimate (2651 in the corresponding proof. The fact that Ux G 



under the same hypotheses follows from the representation (obtained by differentiating ( 229 1 with respect 
to x) 

Ux{x, t) = ^F'{x + t)+ ^-F\x -t) + ^G{x + t)~ ^-G{x - t) 



(266) 



1 
2 

1 
2? 







[b"ix + it-to))-b"ix-it-to))] dto 



[sin(?7(x +{t- to), to) + b{x + (t - to))) 

- sin(C/(x - (i - to), to) + b{x ~{t- to)))] dto , 

which is analyzed precisely in the same way as Ut was in the proof of Theorem |B.5| leading to an estimate 
of exactly the same form as ( 265 ) . 
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Now we suppose further that F" G and G' £ L''{M.), and consider the formula (obtained by 

differentiating (261 1 and (2661 with respect to x: 

^Utjx, t) = ]^F"{x + t)T \f"{x -t) + ]^G'{x + t)± ]^G'{x ~ t) 



(267) 



2^ r(C+)±&"'(e_)]dto 



2 

1 

2? 



[cos(C/(e+, to) + 6(e+)) • (C/.(e+, to) + &'(e+)) 

± cos([/(e-, to) + 6(e-)) • (C/x(e-, io) + &'(?-))] dto , 
where ^± := x ± (i — to)- By Minkowski's inequahty and using | cos(a;)| < 1, 



d_ 
dx 



LP 



\G'\ 



LP 



LP 



(268) 
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1 
2? 

1 
2? 





1 




H 

LP 


" 2 


i 



6"'(- + (t-to))rfto 
|6'(- + (t-to))|rfto 
|f/.(- + (t-to),to)|rfto 



b"'{--{t-to))dto 



LP 





1 




H 

LP 


' 2? 


i 



|6'(--(t-to))|dto 



LP 





1 




H 

LP 


' 2? 


I 



\U,{--{t-to),to)\dto 



LP 



Now, by Holder and Fubini arguments, 



/ |f/,(-±(t-to),to)|dto 
Jo 



LP 



+ C30 



|t/,(a;±(t-to),to)|rfto dx 



(269) 



1 



1 

ti 



+00 pt 



\U,{x±{t~to),to)\Pdto dx 



|l/,(-,to)|lipdto 



<t sup ||f/.(-,to)||LP 
0<to<t 



Applying the same argument to the remaining integrals yields the estimate 
d 

-I I .i- f\ < II ^"'11 r„ -I- lir^'ll r„ -I- / r„ -I I I /l' I I r „ -I- - 

c2 



(270) 



ox X 



< ll^^"llLP + i|G'|Up + <||fe"'||LP + 4ll^'IU- + 4 IIc^.(-.^o)IIlp. 

LP ^ ^ 0<t<to 

Since it has already been established that Ux is in as a function of x uniformly for t in bounded intervals, 
the proof is complete. □ 



The case most relevant for inverse-scattering theory is p = 1. Here we have the following result. 



Corollary B.7. Suppose the initial data for the Cauchy problem (223 1 satisfy sin(/), 1 — cos(/), /', /", g, g' G 
. Then there is a unique global weak solution of the problem (223 1 for which sin(u),l — cos{u),Ux, 
Uxx,ut,utx e for all t>0. 

Proof. We need to translat e the given conditions on / and g into corresponding conditions on F and G 
sufficient to apply Theorem B.5 and Theorem B.6 Since f G L^, f is absolutely continuous and uniformly 
bounded, and therefore so is F. Moreover, the limits F{zLoo) both exist and vanish. The condition that 
sin(/) is in therefore guarantees (since | sin(a; — 27rn)| > \x\/2 for all n G Z and for all a; G M sufficiently 
small) that F G L^. Next, since F'{x) = f'{x) - b'{x) and F"{x) = f"{x) - b"{x), the triangle inequahty 
shows that the conditions f',f" G easily imply that F',F" G L^. Finally, since g differs from G only by 
a factor of e, g,g' G implies G, G' G L^. 

From Theorem B.5 and Theorem B.6 we therefore obtain that U{-,t), Ux{-,t), Uxx{',t), Ut{-,t), and 
Utx{',t) all lie in for all t > 0. Since ut = Ut and Utx — Utx, the fact that Ut, Utx G guarantees that 
ut, Utx G L^. Since Ux — Ux + b'{x) and Uxx = Uxx + b"{x), the triangle inequality shows that Ux, Uxx G 
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guarantees that u^^u^x G L^- Finally, since | sin(w(a;, t)) | = |sin(J7(x,t) + b{x))\, the inequality (2461 shows 
that U E implies that sin(u) £ L^, and the inequality (like (2461, but using instead 1 — cos(a;) < |x|) 

1 — cos(u(a;, t)) — 1 — cos{U{x, t) + b{x)) 

(271) = [1 - CQs{U{x,t) + b{x))]x-{x) + [1 - cos(C/(x,t) - b{~x))]x+{x) 

< \U{x,t)\ + c{x) 

shows that U E implies that 1 — cos(u) G L^, which finishes the proof. □ 

In fact, for the p = 1 case it is possible to show further that an arbitrary number of a:-derivatives 
of u and ut are in L^(M) for all t > if the same holds true at i = 0. The issue in obtaining higher-order 



regularity for general p arises from replacing d/dx with 9" /9a;" in (2671, because while 

d 



(272) 



dx 



sin(/) 



-|cos(/)/'|<|/'| 



gives an estimate that is linear in already-estimated derivatives, the corresponding estimate of the nth-order 
partial derivative will contain, in addition to a term |/^"^ | , a sum of nonlinear terms in lower-order derivatives. 
In the case of p = 1, /('^"^^ is controlled in L°° by by the fundamental theorem of calculus, so all 

of the nonlinear terms may be estimated in by peeling off an appropriate number of uniformly bounded 
factors. For example, to analyze U^xx or Utxx, one replaces d/dx with d^/dx^ in (2671 and then it is required 
to estimate the LP{M.) norm of a term like 

(273) J -^sm{U{x±{t~to)M) + h{x±{t-to)))dh. 

But, since 

• ' = |sin(/)(/')' + cos(/)/"| < iff + |/"| 



(274) 



dx"^ 



sin(/) 



the quadratic term would cause some difficulty for general p. However, for p = 1, the knowledge that 
/" e allows one to further estimate the right-hand side by ii'l/'l -I- |/"| for some constant K that depends 
on Then using /' € as well, the argument proceeds as in the proof of Theorem B.6 

concludes ultimately that Uxxx and Utxx are also in for all i > 0. This general method valid for p 
allows all x-derivatives of Ux and Ut to be handled in the same way. 



and one 
1 
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